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THE  ENERGY  BUDGET  AT  THE  EARTH’S  SURFACE 

Part  II 

INTRODUCTION 


Unfortunately,  the  evaluation  of  the  partition 
of  solar  radiation  into  other  forms  of  energy  at  the 
earth’s  surface  has  to  depend,  in  part,  upon  indirect 
means.  For  example,  energy  exchange  calcula- 
tions based  upon  aerodynamic  principles  have  to 
be  used  to  estimate  sensible  heat  exchange.  The 
principles  are  customarily  used  also  to  estimate 
latent  heat  and  momentum  exchange.  More  re- 
cently, the  photochemical  fixation  of  carbon  di- 
oxide in  the  process  of  photosynthesis  has  been 
evaluated  from  aerodynamic  exchange  calculations 
of  carbon  dioxide  transfer  to  vegetated  surfaces. 

Historically,  the  aerodynamic  exchange  princi- 
ples were  first  developed  in  mnd  tunnel  and  lab- 
oratory fluid  flow  studies,  then  applied  to 
atmospheric  exchange  studies  at  the  earth’s  surface 
either  over  water  or  over  land  where  the  vegetation 
was  generally  short  or  sparse.  It  has  been  com- 
monly accepted  that  such  calculations  are  more 
complicated  over  water  because  of  its  smface  in- 
stability. This  arises  from  the  fact  that  in  the 
exchange  calculations  commonly  used,  it  is  neces- 
sary to  know  certain  boundary  conditions,  i.e., 
have  knowledge  of  a reference  plane  where  the 
turbulent  boundary-  layer  begins.  One  can  readily 
see  that  on  land  surfaces  the  problem  should  be 
much  simpler,  since  the  stable  soil  surface  can  serve 
as  the  reference  plane. 

Recent  investigations  over  more  lush  vegetative 
surfaces  of  the  land,  however,  have  demonstrated 
that  these  surfaces  may  theoretically  and  experi- 
mentally be  more  difficult  to  handle  than  water 
surfaces  simply  because  of  changing  geometric 
characteristics  (plant  growth,  leaf  faU)  and  chang- 
ing elastic  properties  (stalk  bending,  leaf  flutter, 
and  streamlining) . Of  course,  these  problems  are 
amphfied  with  larger  plants;  i.e.,  prairie  grass, 
farm  crops,  and  forests. 

When  one  wants  to  evaluate  the  heat  budget 
over  extended  periods  of  time — for  days  and 
rnonths  as  we  are  attempting  to  do — it  is  prac- 
tically impossible  at  present  to  use  any  aerody- 
namic technique  other  than  the  so-called  log 
profile  methods.  Since  this  method  depends  upon 
knowledge  of  the  surface  reference  plane  where 
the  turbulent  boundary  layer  begins  (and  we  are 
particularly  interested  in  the  energy  exchanges 


with  extensive  vegetative  surfaces),  it  is  necessary 
to  resolve  the  new  complications  that  have  nat- 
urally developed.  This  requires  a systematic 
research  effort  into  the  turbulent  exchange  charac- 
teristics within  and  above  the  dynamic  vegetative 
canopy.  This  report  takes  up  in  detail  the  theo- 
retical and  experimental  approaches  used  in  1960 
to  understand  further  the  aerodynamic  exchange 
processes  within  and  above  the  plant  canopy. 
The  studies  were  carried  out  in  addition  to  the 
routine  collection  of  heat  budget  data  for  the  1960 
growing  season  over  an  alfalfa  field  and  a cornfield. 
The  results  of  the  heat  budget  studies  will  be 
taken  up  in  another  publication. 

The  first  three  sections  in  the  present  study 
take  up  the  theoretical  development  and  support- 
ing experimental  results  concerning  the  airflow 
characteristics  within  the  source  and  sink  region 
of  the  turbulent  boundary  layer  at  the  earth’s 
surface,  as  well  as  the  airflow  characteristics  just 
above  the  source  and  sink  region.  Two  models 
for  “canopy  flow”  (in  the  source  and  sink  layer) 
were  developed.  It  was  later  found  that  the 
model  based  upon  fully  turbulent  flow  within  the 
plant  canopy  most  realistically  applies  to  nature. 
To  facilitate  experimental  studies,  it  was  postu- 
lated that  the  turbulent  structure  could  be  broken 
down  into  two  classifications:  (1)  a transient  state 
and  (2)  a quasi-steady  state.  To  date,  only  the 
quasi-steady  state  has  been  investigated.  The 
turbulent  canopy  flow  model  postulated  a linear 
growth  of  diffusivity  within  the  canopy  flow  layer, 
rather  than  the  log  velocity  profile  (which  holds 
above  the  canopy  flow  layer).  The  theory  and  ex- 
perimental results  demonstrate  a curvature  re- 
versal of  the  velocity  profile  deep  in  the  canopy 
layer. 

The  last  two  sections  discuss  the  complexity  of 
the  relationship  of  the  logarithmic  velocity  profile 
characteristics  above  the  canopy  to  the  dynaniic 
surface  properties  of  lush  vegetation.  To  facili- 
tate the  study  of  the  relationships  involved,  a 
machine  computation  method  was  worked  out  and 
described  (p.28).  From  the  method  one  can  eval- 
uate the  friction  velocity,  roughness  length,  and 
effective  displacement  parameters  along  with  their 
standard  errors.  The  last  section  takes  up  the 
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experimental  findings  from  log  profile  studies  over 
alfalfa,  wheat,  and  corn.  An  amazing  phenom- 
enon was  discovered  in  that  both  the  alfalfa  and 
wheat  demonstrated  a constant  friction  velocity 
with  increasing  windspeed  once  a ‘ ‘critical”  wind- 
speed  was  reached.  Evidently,  the  elastic  plants 
“streamlined”  so  that  a decrease  in  roughness 
compensated  for  increasing  windspeed. 

The  last  section  also  takes  up  an  evaluation  of 


the  distribution  of  shear  and  the  momentum  trans- 
fer coefficient  within  the  vegetative  canopy  in 
wheat  and  corn.  These  studies  focus  attention 
on  the  fact  that  shear  originates  within  the  rough- 
ness volume,  and  not  at  some  empirical  plane. 
Further  studies  are  to  be  made  to  evaluate  also 
the  source  and  sink  distribution  for  sensible  heat, 
latent  heat,  and  carbon  dioxide  within  the  canopy 
layer. 


EFFECTS  OF  TURBULENCE  AND  PHOTOSYNTHESIS  ON  CO 
PROFILES  IN  THE  LOWER  ATMOSPHERE 

By  D.  E.  Ordway,  A.  Bitter,  D.  A.  Spence,  and  H.  S.  Tan 


BACKGROUND 

The  two  causes  of  turbulence  in  the  atmosphere 
are  (1)  forced  convection,  through  the  boundary 
layer  type  shear  flow  instability  (effect  of  velocity 
profile) , and  (2)  free  convection,  through  the  buoy- 
ancy force-type  instability  (effect  of  heating  from 
below) . 

At  present  much  literature  on  the  study  of 
turbulence  exists  on  forced  convection,  but  little 
deflnitive  knowledge  is  available  regarding  the 
second  type.  In  the  present  study,  it  is  our  pur- 
pose to  limit  the  first  phase  of  the  investigation 
to  the  effect  of  the  forced  convective,  or  boundary- 
layer  type  turbulence  only,  on  atmospheric  CO2 
distribution.  We  are  considering  the  problem  of 
steady  state  atmospheric  boundary  layer  flow  and 
shaU  now  introduce  a principle  of  “flux  conserva- 
tion.” 

This  principle  states  in  essence  that  the  trans- 
port equation  for  any  physical  quantity  Q must 
be  of  the  form 

^=^+2.ve=div(iicve)+sso  (1) 

where  K is  the  flux  coefficient,  S the  distributed 
source  distribution,  t is  tune,  and  g the  flow 
velocity.  In  other  words,  particle  derivative  of 
Q must  vanish.  This  is  indeed  a generalization  of 
the  well-known  constant  flux  hypothesis  {29)  ^ 
long  adopted  by  micrometeorologists,  to  include 
the  distributed  source  contribution. 

Under  the  usual  hypothesis  of  flow  stratifica- 
tion, application  of  this  principle  to  the  steady 
state  distribution  of  CO2  concentration,  E,  with 
height,  z,  immediately  gives 

2.  PV£=|  [/>!?] (2) 

where  w is  distributed  sink  or  source  of  CO2  due 
to  photosynthesis  and  respiration,  and  D the  total 
diffusivity;  i.e.,  sum  of  molecular  diffusion  (lam- 
inar, Dl)  and  eddy  mixing  (turbulent,  Dt)  . 

Application  of  the  same  principle  to  x-momen- 
tum  with  u as  wind  velocity  gives 

?.  (3) 


where  ^ is  distributed  shear  source,  i.e.,  leaf  drag, 
and  T is  the  shear  stress,  given  by 

r = (p  + .)g  (4) 

with  ^ denoting  laminar  viscosity  and  €,  eddy 
viscosity. 

The  validity  of  this  principle  can  most  easily 
be  seen  from  the  limiting  behavior  of  equation  3. 
Indeed,  for  outer  flow  where  t?  = 0,  this  equation 
leads  to  a constant  shear  requirement  that  is 
famous  for  the  logarithmic  velocity  profile  in 
turbulent  wall  flow  and  linear  velocity  profile  in 
laminar  Couette  flow.  On  the  other  hand,  when 
^7^0,  this  equation  directly  gives  the  differential 
equation  of  canopy  flow  velocity  distribution. 

Equations  2 and  3 will  form  the  basis  of  our 
mathematical  study  of  the  atmospheric  carbon 
dioxide  concentration  and  the  canopy  flow,  respec- 
tively. 

ATMOSPHERIC  CARBON  DIOXIDE 
CONTENT  STUDY 

CO2  content  in  the  lower  atmosphere  controls 
and  is  in  turn  controlled  by  the  life  process  on 
earth.  The  mechanism  that  keeps  local  CO2  con- 
tent from  becomiug  deficient  or  excessive  is  the 
process  of  molecular  (laminar)  and  eddy  (turbu- 
lent) diffusion.  The  differential  equation  gov- 
erning the  steady  state  CO2  concentration,  E,  as- 
suming adequate  CO2  for  saturation  requirement 
of  photosynthesis,  can  therefore  be  written  in 
general  as  follows: 


(5) 

w= 

---j{z)E+g{z) 

(6) 

where  D denotes  the  total  diffusivity,  w the  rate  of 
CO2  assimilation  or  production  as  a result  of  photo- 
synthesis (J)  and  respiration  {g).  Functions/  and 
g are  defined  by: 

/(s)  = photosynthetic  rate  coefficient=UiA(2)cr'(2) 

(7) 

1 Italic  figures  in  parentheses  refer  to  Literature  Cited, 
p.  48. 
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^(2)=rate  of  respiration =C20-' (3).  (8) 

Here  L{z)  is  the  light  intensity  function  (5),  (j'{z) 
the  leaf  area  density  function,  and  Ci  and  C2  are 
constants. 

With  w=^  for  z^h  (^=plant  height),  it  is 
evident  that  the  plane  z=h  separates  the  atmos- 
phere into  two  regions,  inside  each  of  which  dif- 
ferent laws  govern  the  CO2  distribution,  namely. 

The  outer  flow,  — • 


The  inner  flow,  z<Ji — 

Equations  5 and  6 show  clearly  the  respective 
roles  played  by  turbulence  and  photosynthesis  on 
the  atmospheric  CO2  distribution:  the  effect 
of  turbulence  appears  in  the  diffusive  term 

through  D,  the  eddy  diffusivity  coeffi- 
cient. The  effect  of  photosynthesis  that  depends 
on  leaf  area,  light  intensity,  and  atmospheric  CO2 
concentration  appears  as  the  linear  term  j{z)E  of 
the  differential  equation.  The  nonhomogeneous 
term  g{z)  represents  the  effect  of  respiration,  which 
is  independent  of  light  and  CO2  contents. 

Now,  the  kinetic  transport  mechanism  for  mass 
(diffusion),  momentum  (shear),  and  energy  (heat 
conduction)  are  identical.  Upon  introducing  the 
concept  of  constant  turbulent  Schmidt  number 
defined  by 


(11) 


we  can  usually  establish  some  correlation  between 
CO2  concentration  profile  E{z)  and  velocity 
profile  u{z).  Let  us  consider  now  the  two  flow 
regimes. 


However,  one  must  be  careful  in  not  pushing 
the  analogy  too  far.  Equation  12  is  indeed 
identical  with  that  governing  the  velocity  dis- 
tribution, but  the  analogy  stops  here.  Whereas 
the  nonslip  condition  over  the  datum  plane  of 
velocity  profile  necessarily  represents  a momentum 
sink,  the  boundary  condition  over  the  datum 
plane  of  concentration  profile  can  be  either  source 
or  sink,  depending  solely  on  the  relative  rate  of 
respiration  and  photosynthesis.  Moreover,  even 
the  datum  planes  for  the  two  cases  are  usually  at 
different  locations.  Thus,  the  only  conclusion 
one  can  safely  draw  from  the  similitude  of  the  two 
differential  equations  is  that  the  two  solutions 
must  belong  to  the  same  family  of  curves. 

Inner  Flow, 

In  this  region,  r 5^  constant.  As  can  be  seen  by 
comparing  with  the  canopy  flow  equation  derived 
in  the  discussion  of  the  laminar  flow  model  study, 
p.  5,  the  differential  equation  for  E differs 
essentially  from  that  for  u.  Thus  little  correlation 
between  the  two  profiles  can  be  expected.  Here, 
two  cases  should  be  considered  separately,  de- 
pending on  the  linear  term  in  the  differential 
equation. 

1.  Respiration  Profile  {Nighttime). 

In  this  case, 

j{z)  =L{z)  =0,  g{z)  5^0.  (14) 

Equation  10  reduces  to  the  following  form, 

A(z,g)+,(.)=0.  (15) 


The  solution  to  this  equation  can  be  obtained 
through  integration,  as  follows: 


dz 

E-Eo= 


(16) 


Outer  Flow,  z)>h 

Since  outer  flow  is  nothing  but  the  turbulent 
wall  flow,  where  t=to,  e—Cz,  on  introducing  the 
Schmidt  number,  equation  9 takes  the  form 


giving 


dzV  bz)  ° 

(12) 

E—EQ=C\n.  z. 

(13) 

Like  the  velocity  distribution,  E also  has  a 
logarithmic  profile. 


2.  Photosynthesis  Profile  {Daytime). 

In  this  case,  if 

J{z)9^0,  g{z)9^0, 

the  full  equation  10  must  be  retained;  i.e.. 


^ d^E  dD  dE  r/  X 77  / N 


(17) 

(18) 


A general  solution  to  this  equation  usually  pre- 
sents great  complication  and  difficulty.  How- 
ever, if  we  make  the  approximation  P= constant, 
which  has  been  commonly  adopted  by  various 
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investigators  {32),  this  equation  simplifies  to  the 
following  generalized  Hille’s  equation; 

D^-J{B)E=-g{z).  (19) 

Although  this  assumption  is  admittedly  crude, 
nevertheless,  it  should  provide  a good  first  ap- 
proximation. 

A solution  to  equation  19  can  be  obtained  by 
introducing  a Fourier  spectrum  for  J{z)  and  ex- 
pressing the  complementary  solution  in  the  form 

(20) 

The  determination  of  characteristic  exponents 
by  Hille  determinant  and  the  recursion  formulas 
connecting  the  a^’s  and  bjs  are  discussed  in  detail 
by  Whittaker  and  Watson  (39).  For  the  partic- 
ular solution,  by  witing 

Y^v,Y^=v,Y,  (21) 

and  applying  the  method  of  variation  of  para- 
meters, it  can  be  shown  that  (10) 

where  A(Fi,F2)  is  the  Wronskian  of  Yi  and  F2. 

CANOPY  FLOW  VELOCITY 
DISTRIBUTION  STUDY 

We  define  canopy  fiow  as  the  flow  inside  the 
plants  {0<C.z<^h),  which  consists  mainly  of  a 
turbulent  layer  at  the  upper  part  and  a wake-type 
fiow  at  the  lower  part.  It  is  clear  that  a really 
satisfactory  study  must  be  a statistical  one,  of 
ensemble  average  type  approach.  However,  as  a 
crude  approximation,  it  seems  also  plausible  that 
we  can  replace  the  distributed  plant  leaves  by  a 
fictitious  distribution  of  shear  centers  and  solve 
the  differential  equation  corresponding  to  this 
model.  This  approach  is  by  nature  speculative. 
At  the  present  moment  there  are  two  models  of 
approach  under  consideration. 

Laminar  Flow  Model 

The  differential  equation  for  velocity can  be 
written: 

flu'll 

g~F{z)u^=0  (23) 


where  is  the  laminar  viscosity  and  F{z)  is  some 
function  of  the  plant  characteristics.  The  indi- 
rect solution  to  this  equation  can  be  obtained  by 
quadrature  after  certain  transformations  and  ex- 
pressed in  terms  of  elliptic  functions  as  follows: 

(^J=?  £ F[ziu)]u^du+kt 

SO, 

A detailed  discussion  and  some  typical  results  of 
this  model  are  given  below. 

In  the  first  model  of  the  canopy  flow  velocity 
distributions  study,  an  approximate  solution  to 
the  velocity  profile  u{z)  in  the  base  flow  within  the 
plants  has  been  considered.  At  the  upper  edge 
of  the  flow,  say  z = zi,  we  see  that  the  outside 
turbulent  flow  effectively  “drags”  the  canopy 
flow,  analogous  to  a moving  wall  with  a velocity 
Ui.  If  we  replace  the  leaves  by  an  effective  distri- 
bution of  infinitesimal  “shear  centers,”  we  can 
write  an  equation  of  horizontal  force  equilibrium 
from  figure  1.  We  assume  that  r,  the  fluid  dy- 
namic shear,  is  due  mainly  to  the  mokcular 
momentum  transport  of  laminar  flow  and  D,  the 
average  drag  per  unit  volume  of  the  shear  centers, 
is  turbulent  in  nature.  In  nondimensional  form, 
this  yields, 

0=g_S(5)  I (25) 

where  u = uju  i,z  = zlzi  and 

S(z)  = CoRe,'^^-  (26) 

Cd  is  a mean  drag  coefficient  of  the  shear  centers, 
which  is  turbulent  and  is  independent  of  u;  Rei  is 
the  Reynolds  number  dependent  on  Ui  and  Zi,  and 
da(z)ldz  is  an  effective,  equivalent  local  leaf  area 
density.  Equation  25,  together  with  the  boundary 
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conditions  at  Zi  and  0,  respectively, 

u(l)  = l (27) 

u{0)=0  (28) 


Further  study  of  both  models  should  be  carried 
out.  A final  choice  can  only  be  made  after 
sufficient  field  data  have  been  gathered  to  confirm 
the  theoretical  model. 


then  determines  the  velocity  profile  in  the  canopy 
flow.  Since  the  solution  is  characterized  only  by 
S{T),  we  call  S{z)  the  “shape  function.”  Because 
it  is  always  positive,  we  see  that  the  variation  of 
the  slope  is  always  negative.  That  is,  the  velocity 
profile  is  always  concave  backwards  as  opposed 
to  the  typical  laminar  boundary  layer,  which  is 
concaved  forward. 

Equation  25  is  a quasi-nonlinear  differential 
equation  of  second  order.  In  general,  the  solu- 
tion has  not  been  found.  However,  for  constant 
S or  S{z)=S,  a suitable  transformation  to  the 
velocity  gradient  as  the  independent  variable  and 
the  velocity  as  the  dependent  variable  permits 
integration.  The  result  is 


z=F  \ cos 


{ylZ-l)k-U 

( V3  ~h  1 ) ^ “1“ 


2+V3J 


-^’[cos->  (2-V3),  (29) 


where  k is  related  to  S and  the  profile  slope  at  the 
ground;  i.e.,  the  ground  shear  by 


and  F is  the  elliptic  function  of  the  first  kind. 
For  a given  value  of  the  constant  shape  function, 
the  parameter  k,  and  hence  the  ground  shear,  is 
determined  in  general  by  satisfying  equation  27. 
The  results  are  plotted  in  figure  2 over  a repre- 
^ntative  range.  A typical  velocity  profile  for 
^=141  has  also  been  computed  and  is  illustrated 
in  figure  3. 

Turbulent  Flow  Model 


This  model  has  been  suggested  under  the  heu- 
ristic argument  that  on  the  leaf  scale,  the  fiow  may 
appear  laminar;  on  plant  scale,  it  must  appear 
turbulent.  The  corresponding  differential  equa- 
tion takes  the  form: 


where  e{z)  is  the  eddy  viscosity. 

For  a first  approximation,  we  may  take  e as 
constant  (several  orders  of  magnitude  greater 
than  jLt  (55)),  thus  reducing  equation  31  to  equation 
24,  only  with  a different  scale. 


k 


Figure  2. — Canopy  flow  wall  shear  for  a constant  shape 
function. 


Relative  Height  (z) 


Figure  3. — Canopy  flow  velocity  profile. 


QUASI-STEADY  MICRO-METEOROLOGICAL  ATMOSPHERIC 
BOUNDARY  LAYER  OVER  A WHEATFIELD 

By  H.  S.  Tan  and  S.  C.  Ling 


THEORETICAL  MODEL 

Concept  of  Relevant  Scale 

In  the  study  of  turbulent  transport  process  in 
the  lowest  part  of  the  earth’s  atmosphere,  one  of 
the  most  important  concepts  is  the  “relevant 
scale”  of  the  micro-structure  of  the  wind.  Indeed, 
any  typical  wind  record  taken  at  a point  in  a corn- 
field or  wheatfield  exhibits  (fig.  4)  at  least  two  dis- 
tinct ranges  of  eddy  scales:  (1)  a high  frequency 
fluctuating  velocity  component  of  small-scale 
eddies  superimposed  on  (2)  a lower  frequency 
gusty  wind  composed  of  large  eddies.  To  a cosmic 
observer  in  the  large  or  meteorological  sense,  both 
these  eddies,  with  scales  of  different  orders  of  mag- 
nitude, are  nothing  but  turbulence.  However,  to 
an  observer  sitting  at  the  top  of  a wheat  crop  in 
the  small  or  micro-meteorological  sense,  the  pic- 
ture is  different.  What  he  recognizes  here  is  a 
turbulent  (high  frequency),  gusty  (low  frequency) 
wind.  The  recognition  of  such  a “relevant  scale” 
is  important  if  we  recall  that  the  term  “turbulence” 
is  so  vaguely  defined.  Indeed,  P.  A.  Sheppard 
brought  this  point  out  most  clearly  in  his  final 
address  to  the  International  Symposium  on  At- 
mospheric Turbulence  in  the  Boundary  Layer  (7) . 
He  points  out  as  follows : 

The  atmosphere  . . . shows  simultaneously  many 
modes  of  motion,  varying  in  scale  from  the  “long 
waves”  of  the  westerlies  down  through  depressions 
and  anticyclones,  convective  motions  and  so  on,  to 
small  scale  turbulence  . . . 

Now  the  motions  which  we  commonly  study  under 
the  name  of  turbulence  are  so  much  smaller  in  scale  than 
depressions  and  the  like  that  we  make  little  or  no 
attempt  to  pick  out  individual  elements  of  the  flow, 
except  quite  formally  as;  e.g.,  in  a power  spectrum. 
But  who  doubts  that  if  we  magnified  the  fleld  suffi- 
ciently, drew  a series  of  “synoptic  charts,”  in  fact, 
we  should  in  general  find  recognizable  structures  and 
changes  in  those  structures  with  time  . . . 

If  we  study  the  wind  structure  in  its  spectrum 
over  a relatively  short  time,  we  notice  that  the 
energy  of  turbulence  usually  resolves  itself  into 
two  widely  separated  wave  packets  (fig.  5).  This 
phenomenon  can  be  interpreted  from  a micro-mete- 
orological point  of  view  as  follows:  Ei,  associated 
with  small-scale  eddies,  carries  the  turbulent 
mixing  part  of  the  kinetic  energy  that  controls 
both  the  vertical  and  horizontal  transport  proc- 
esses; E2  associated  with  large-scale  eddies. 


0 50  100  150 

Time  (sec.) 


Figure  4. — Horizontal  wind  (u)  at  7-meter  height. 


Figure  5. — Typical  spectrum  of  horizontal  wind. 

represents  the  part  of  the  kinetic  energy  carried 
away  by  gusty  winds,  which  can  only  affect  the 
horizontal  transport.  This  separation  of  two  dis- 
tinct eddy  scales  responsible  for  the  vertical  and 


8 


PRODUCTION  RESEARCH  RPT.  7 2,  U.S.  DEPT.  OF  AGRICULTURE 


the  horizontal  transport  processes,  respectively.  Now  let  X denote  the  eddy  size,  or  wavelength, 
was  also  arrived  at  by  A.  Monin  (20).  and  then  we  have  the  following  scheme: 


Subdivision  of  Atmospheric  Boundary  Layer : 

Synoptic  Versus  Micro -Layers 

Whenever  there  is  a fluid  flow  past  a surface,  a 
boundary  layer  will  develop.  Now  the  structure 
of  the  atmospheric  wind  consists  of  components 
of  large-  and  small-scale  eddies  embedded  in  the 
synoptic  (geostrophic)  wind.  The  large-scale 
turbulence  present  in  the  synoptic  wind  boundary 
layer  is  caused  by  the  presence  of  large-scale 
roughness  elements,  such  as  terrain  irregularities, 
and  the  layer  thickness  usually  extends  to  a height 
of  a few  kilometers.  The  relatively  smaller  scale 
micro-meteorological,  or  gusty  winds,  are  nothing 
but  these  large-size  turbulent  eddies  in  the  synoptic 
layer.  Indeed  these  are  actually  air  lumps  of  the 
size  of  a few  hundred  meters,  sweeping  across  the 
field  and  leaving  behind  the  gust  front  a growing 
micro-boundary  layer.  This  layer  quickly  reaches 
a quasi-steady  state  a few  meters  behind  the  prop- 
agating front  and  resembles  very  much  the  steady 
flow  leading  edge  boundary  layer  over  a flat  plate. 
The  thickness  of  this  micro-boundary  layer  inside 
the  micro-meteorological  wind  usually  extends  to 
about  10  meters.  Further,  the  small-scale  tur- 
bulence inside  this  layer  with  eddy  sizes  ranging 
in  decameters  is  caused  by  small-scale  roughness 
elements,  such  a plant  crops. 

On  the  average,  a gusty  wind  of  10  meters  per 
second  with  a steady  interval  of  several  minutes 
and  a turbulent  fluctuating  component  having  a 
frequency  of  a fraction  of  a cycle  per  second  are 
typical.  On  the  other  hand,  a micro-wind  in  a 
field  is  never  steady.  However,  the  relatively 
well-defined  edge  of  gusts  (air  lumps)  indicates 
that  a frontal  theory  can  be  developed.  The 
highly  transient  portion  behind  the  gust  front 
extends  about  10  meters,  after  which  a quasi- 
steady state  is  quickly  established.  This  means 
that  a 10-second  average  field  measurement  taken 
in  the  middle  portion  of  a 1 -minute  period  gusty 
wind  should  indicate  a consistent  steady  state 
logarithmic  velocity  profile  outside  the  crop. 
The  validity  of  this  scaling  idea  has  now  been 
amply  verified  by  field  measurements. 

The  classification  of  eddy  sizes  and  definition  of 
a quasi-steady  condition  can  be  quantitatively 
expressed  in  terms  of  the  micro-boundary  layer 
thickness,  8,  and  the  relaxation  time,  t,  in  the 
following  manner. 

Quasi-steady  condition  considered  reached  at: 


Scale 

^<2,0.... 

2.0<^<20 

^>20..... 


Eddy  classification 
Turbulent  eddies. 

Unsteady  gusts. 

Quasi-steady  wind. 


Proper  Averaging  Period 

In  a natural  field,  over  a sufficiently  long  period, 
all  scales  of  eddies  will  appear  and  the  spectrum 
will  be  continuous.  A velocity  mean  taken  over 
an  arbitrary  period  of  time  is  thus  meaningless. 
The  proper  averaging  period  must  be  arrived  at 
by  “relevant  scale”  considerations  in  accordance 
with  the  process  involved.  For  example,  in  the 
vertical  transport  process  that  is  controlled  solely 
by  small-scale  turbulence  (wavelength  of  meter 
size  and  period  of  seconds),  extraneous  effects  of 
the  large-scale  turbulence  (gusty  winds)  will  be 
introduced  if  the  averaging  period  is  taken  in 
hours. 

This  will  utterly  ruin  the  statistical  behavior  of 
the  relevant  small-scale  eddies.  On  the  other 
hand,  for  the  horizontal  transport  process  that  is 
governed  by  all  scales  of  turbulence,  averaging 
periods  of  less  than  a minute  will  give  no  trend  at 
all.  Moreover,  the  results  will  be  meaningless 
since,  under  such  an  averaging  process,  the  largest 
energy-carrying  eddy  scales  in  the  spectrum  are 
completely  left  out. 

The  technique  proposed  in  the  present  report 
therefore  involves  a selective  sampling  procedure, 
in  order  to  insure  the  proper  quasi-steady  be- 
havior. It  must  be  remarked  at  this  point  that 
this  quasi-steady  state  study  alone  is  not  sufficient 
for  complete  knowledge  of  the  vertical  transport 
process;  this  must  be  augmented  by  a further  study 
of  the  transient;  i.e.,  an  accelerating  and  decelerat- 
ing regime.  The  final  effective  eddy  diffusivity 
will  have  to  be  evaluated  as  a weighted  sum  of  the 
quasi-steady  and  nonsteady  measurements,  and 
should  more  or  less  deviate  from  that  correspond- 
ing to  the  classical  logarithmic  velocity  profile. 
For  the  study  of  the  nonsteady  regime,  frontal 
■heory  might  be  expected  to  play  an  essential 
role.  However,  such  a study  will  not  be  attempted 
before  we  enter  the  second  phase  of  our  research 
program. 


?7o^>5  where  5=  10  m.. 


FIELD  MEASUREMENTS 


^>yySeC. 

^ = 0 at  the  gust  front. 


Experimental  Considerations 

As  stated  earlier,  in  a natural  field  over  a 
sufficiently  long  time,  the  entire  spectrum  of  eddy 
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scales  appears  and  hence,  theoretically,  the 
relevant  scale  encompasses  the  complete  spectrum. 
Such  longtime  averages  of  atmospheric  measure- 
ments are  definitely  of  interest  in  the  study  of 
weekly  or  seasonal  variations  of  atmospheric 
phenomena  and,  in  particular,  of  horizontal 
transport  processes.  On  the  other  hand,  over  a 
relatively  short  interval,  we  can  usually  identify 
two  distinct,  dominating  ranges  of  eddy  scales,  of 
which  the  smaller  one  alone  controls  the  vertical 
transport  process  in  the  earth’s  atmosphere. 
Thus,  inasmuch  as  we  restrict  our  study  to  the 
vertical  atmospheric  distribution,  the  long  period 
averages  can  be  ruled  out  as  improper.  The 
‘‘relevant  scale”  for  the  vertical  distribution 
study  is  the  small  one;  the  corresponding  averaging 
time  is  the  short  periods;  i.e.,  of  the  size  of  meters 
and  period  of  10  seconds. 

Based  on  the  above  reasoning,  a new  approach 
for  handling  this  problem  was  proposed:  separa- 
tion of  the  transient  part  of  the  wind  from  the 
steady  part  and  a study  of  their  characteristics 
as  separate  entities.  In  order  to  envision  more 
clearly  the  picture  involved,  one  may  consider 
the  atmospheric  boundary  layer  as  consisting  of 
three  distinct  cascading  layers,  as  follows: 

(1)  Laminar  sublayers,  at  the  surface  of  the  leaf 
and  stalk; 

(2)  Micro-meteorologic,  turbulent  boundary 
layer  (over  plant  top  and  canopy  flow  within 
plants) ; 

(3)  Synoptic,  turbulent  boundary  layer  of 
earth’s  atmosphere. 

These  three  layers  embrace  different  scales  of 
eddies  and  have  different  mean  velocity  gradients. 
The  roughness  of  leaf  and  stalk  surfaces  as  seen  by 
the  laminar  sublayer  does  not  directly  influence 
the  mixing  mechanism  of  the  micro-turbulent 
layer.  On  the  other  hand,  the  fluctuating  velocity 
in  the  laminar  sublayer  is  directly  the  result  of  the 
eddies  of  the  micro-turbulent  layer.  In  the  same 
manner  the  roughness  of  plants  as  seen  by  the 
micro-turbulent  layer  plays  no  direct  role  in  the 
mixing  mechanism  of  the  larger  synoptic  atmos- 
pheric boundary  layer.  However,  the  large  eddies 
of  the  synoptic  layer  are  responsible  for  the  gusti- 
ness of  the  micro-turbulent  layer.  This  is  pre- 
cisely why  the  micro-turbulent  layer  is  never 
steady.  The  turbulence  in  the  synoptic  layer  is 
controlled  essentially  by  the  roughness  scale  of 
much  larger  magnitude,  such  as  terrestrial  irregu- 
larities (hflls  and  valleys),  other  disturbances  gen- 
erated by  temperature  gradients  (buoyant  mo- 
tion), and  continental  pressure  fields. 

The  small-scale  eddies  in  the  micro-turbulent 
layer  are  now  solely  responsible  for  the  vertical 
diffusion  of  carbon  dioxide,  vapor,  and  heat  around 
plants.  This  is  the  layer  that  is  to  be  studied  in 
detail,  in  order  to  understand  the  mechanism  of 
lower  atmospheric  diffusion.  As  mentioned  ear- 
lier, this  layer  is  influenced  by  the  large  eddies  of 


the  outer  synoptic  boundary  layer,  and  therefore 
the  state  of  flow  varies  constantly  between  the 
transient  and  the  quasi-steady.-  During  the  tran- 
sient stage,  the  wind  shear  is  not  constant  over  the 
height,  so  that  the  logarithmic  wind  profile  does 
not  exist.  However,  with  eddies  of  sufficient  size, 
a period  of  quasi-steady  state  exists  and  this  then 
enables  one  to  measure  the  steady-state  flow- 
parameters. 

As  a first  step  toward  the  understanding  of  a 
complex  problem,  this  report  is  limited  to  the  study 
of  the  quasi-steady  state. 

The  wheatfield  selected  for  the  present  study  is 
situated  on  top  of  Mount  Pleasant,  5 miles  east  of 
Ithaca,  N.Y.  Measurements  of  wind  profiles 
were  taken  during  July  1960,  starting  when  the 
heads  of  wheat  were  level  with  the  top  leaves,  and 
continuing  up  to  the  time  when  the  heads  were 
30  cm.  above  the  top  leaves.  The  average  peak- 
plant-height  ranged  from  95  to  130  cm.  Average 
plant  density  was  60  stalks  per  square  foot. 

The  instrument  mast  was  located  near  the  center 
of  a 200-  by  200-meter  field.  Small-cup  anemome- 
ters, 9 cm.  in  diameter,  were  used.  They  were 
seated  at  10,  20,  40,  80,  and  160  cm.  above  the 
peak  level  of  crop.  Counters  for  all  anemometers 
were  read  simultaneously  at  every  30-second  in- 
terval. Attempts  were  made  to  select  the  wind- 
profile  record  only  when  a quasi-steady  wind  ex- 
isted. That  is,  a wind  record  was  selected  when 
it  was  preceded  and  followed  by  a 30-second  wind 
record  with  approximately  the  same  average  wind 
velocity  as  measured  by  the  top  anemometer. 

From  the  experimental  record  it  was  evident 
that  for  a considerable  part  of  the  time  the  wind 
profiles  over  the  crop  were  transient. 

Five  representative  steady  wind-profiles  of  dif- 
ferent mean-velocity  ranges  are  plotted  in  figure  6. 
One  important  characteristic  is  that  the  mean  ve- 
locity gradient  at  the  peak  level  of  the  crop  changes 
only  slightly  with  increase  in  mean  velocity.  In 
addition,  the  characteristics  of  these  profiles  are 
insensitive  to  thermal,  conditions  and  stages  of 
plant  growth. 

Logarithmic  Wind  Profiles 

Steady  wind  velocity  profiles  over  a crop  can  be 
expressed  by 

V(Z,  log 

where 

k=0A2,  von  Karman’s  constant; 

Z= height  above  ground; 

A, = peak  level  of  crop; 

t7o=mean  velocity  outside  of  boundary  layer 
(at  10-meter  height  based  on  logarith- 
mic profile) ; 
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U^=ln  (10)- 


, friction  velocity; 


(f= zero-plane  displacement; 
Zo= roughness  length; 

- T= shear  stress; 
p=mass  density. 


, d,  and  Zq  are  to  be  determined  from  the  meas- 
ured wind  profiles,  and  they  are  found  to  be  de- 
pendent on  Uq.  The  value  of  d is  obtained  from 
equation  33,  where  the  mean  velocity  gradients  at 
two  different  heights  are  obtained  from  measure- 
ments. 


is  computed,  after  d is  known,  as  follows: 


u,-u. 


log 


{Z,-d) 

{Z2-d) 


(34) 


Z-h 


Zq  is  computed  after  d and  are  found. 


y Z\  d 

^0— lowi/v/ 


(35) 


When  the  three  parameters  have  been  computed 
for  each  experimental  wind  profile,  the  value  of  JJq 
(at  10  meters  height)  can  be  computed  from 
equation  32. 


Values  of  Zq,  and  ^ are  plotted  in  figure  7 as 


a function  of  JJq.  Both  ^ and  JJ^  are  found  to  be 


insensitive  to  free-stream  velocity.  On  the  other 
hand,  Zq  decreases  rapidly  with  increase  of  JJq  and 
remains  constant  when  JJq  increases  beyond  8 
meters/sec.  The  drag  coefficient 


Cc  [ln(10)P^2  ^(jJ 
2 


is  plotted  as  a function  of  J/q  in  figure  7.  The 
reduction  of  Cz,  with  increasing  J/q  is  evident. 

Canopy  Flow 

The  basic  differential  equation  for  canopy  flow 
(SS)  is 

^h(Z)+,]g=A(Z)^  (36) 

where 


€(Z)=eddy  viscosity,  a function  of  Z; 

/x= dynamic  viscosity,  a constant; 

A(Z)  = total  vertical  plane  area  of  plants/unit 
volume,  a function  of  Z; 

(7^— drag  coefficient  of  plant; 

p=mass  density  of  air; 

77;,= mean  velocity  at  height,  h. 

Within  the  flow  range  considered,  Cd  can  be 
taken  as  independent  of  Reynolds  number. 
Equation  36  can  be  rearranged  as 


— — ULJ 

6(Z)+M  “ 2 


(37) 


where  JJ  and  Z are  JJ  and  Z normalized  by  JJn 
and  h,  respectively.  We  consider  S{Z)  to  be  the 
parameter  (shape  function)  of  the  differential 
equation.  For  wheat,  A{Z)  can  be  considered  as 
independent  of  (Z),  because  the  area  density  of  the 
plant  is  practically  constant  with  height.  The 
eddy  viscosity  e{Z)  may  be  taken  as  being  a 
constant  for  the  top  quarter  of  Z,  for  the  reason 
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that  the  energy-carrying  eddies  are  not  expected 
to  be  modified  extensively  as  they  penetrate  into 
the  plants.  Hence  SiZ)  may  be  taken  as  a con- 
stant S. 

With  boundary  conditions 
^Iz=i=l 


U\z==Q—0 


equation  37  is  unique  and  can  be  solved  for  each 
value  of  S.  The  canopy  wind  profiles  are  plotted 
in  figure  8 for  various  shape  functions  frorn^lO  to 


1,000. 


The  velocity  gradient  at 


bZ 


is 


plotted  in  figure  9. 

The  velocity  gradient  of  the  wind  profiles  at  h 
must  match  the  velocity  gradient  at  the  edge  of 
the  canopy  flow  (fig.  9).  From  this  condition  the 
shape  function,  S,  as  a function  of  outer  boundary- 
layer  velocity,  Uo,  can  be  found  and  is  plotted  in 
figure  7.  This  completes  the  tie  between  the 
wind  profile  and  canopy  flow  once  the  steady  wind, 
Uo,  is  given. 

The  shape  function,  *8,  drops  rapidly  with  in- 
crease of  IJo  and  tends  to  reach  a terminal  value 
above  a 7-meter/sec.  windspeed.  In  figure  8,  one 
notes_that  the  canopy  profiles  are  more  curved  for 
high  S and  low  Uo,  and  approach  asymptotically  a 


Relative  Height  (Z) 


Figure  8. — Canopy  wind  profiles  for  various  constant 
shape  functions,  S. 


1 10  _ 100  1,000 

s 

Figure  9. — Velocity  gradients  at  peak-plant-height  as 
a function  of  shape  function,  S. 


straight  line  for  low  S and  large  Uo.  S varies  di- 
rectly as  CU , and  inversely  as  e(Z);  the  reduction 
of  S at  large  Uo  means  a reduction  in  Cd  and  an 
increase  in  e{Z).  This  in  turn  means  that  the 
flow  is  more  turbulent.  We  see  therefore  that 
the  field  measurements  exhibit  trends  similar  to 
those  anticipated  from  the  theory. 

Accurate  measurement  of  canopy  flow  has  not 
yet  been  made,  but  from  some  preliminary  meas- 
urements there  seems  to  be  a tendency  ^or  the 
canopy  wind  profile  to  cross  from  a lower  S value 
to  a higher  S value  curve,  as  Z goes  below  0.8. 
This  would  mean  that  e(Z^may  not  be  a constant, 
but  should  be  reduced  as  Z is  lowered.  Thus  far, 
experimental  data  are  not  sufficient  to  evaluate 
6(Z). 

From  the  known  value  of  the  S factor,  for  a 
given  Un  and  plant  density,  the  value  of  eddy 
viscosity  can  be  computed.  For  the  range  of  data 
presented,  the  ratio  of  eddy  viscosity  to  dynamic 
viscosity  ranges  from  1 X 10^  to  3X10^.  The  eddy 
viscosity  at  peak  plant  height,  en,  can  be  com- 
puted from  the  logarithmic  equation,  equation  32, 

en  = pk{h-d)U^I\n  (10).  (38) 

The  ratio  of  eddy  viscosity  to  dynamic  viscosity 
obtained  by  this  method  is  considered  to  be  more 
accurate,  and  it  was  found  to  range  from  0.9X10^ 
to  2.2X10^  The  values  at  high  range  do  not 
agree  well  with  those  computed  by  canopy  flow. 
The  error  is  probably  due  to  the  assumption  of 
constant  e{Z).  More  experimental  measurements 
and  analytical  study  on  e{Z)  will  be  conducted. 


A STUDY  OF  ATMOSPHERIC  TURBULENCE  AND  CANOPY  FLOW 


By  H.  S.  Tan  and  S.  C.  Ling 


In  the  present  study  the  following  two  aspects 
are  stressed:  (1)  To  furnish  the  “relevant  scale” 
concept  with  a more  rigorous  foundation  and  (2) 
to  improve  the  canopy"  flow  model  by  postulating 
the  linear  growth  of  diffusivity  rather  than  the 
logarithmic  velocity  profile  as  the  basic  property 
of  turbulent  wall  flow.  It  is  gratifying  to  find 
that  the  predicted  velocity  profile  curvature  re- 
versal deep  in  the  canopy  flow  was  indeed  ob- 
served and  pointed  out  repeatedly  by  several 
investigators  (22). 

It  might  not  be  out  of  place  to  mention  at  this 
point  that  until  now  we  have  not  approached  the 
diffusion  problem  with  source  distribution,  the 
photos3uithesis  problem,  and  the  thermo-insta- 
bility problem  as  seriousty  as  we  would  like  to. 
The  main  reason  behind  this  is  the  heavy  de- 
pendence of  these  phenomena  on  the  atmospheric 
turbulence,  of  which  little  is  definitety  known  and 
much  has  yet  to  be  clarified  before  we  can  satis- 
factority  formulate  a working  model  of  atmospheric 
transport. 

An  introductory’^  discussion  of  turbulent  dif- 
fusivity in  shear  flow  is  provided  by  Dr.  Spence. 

Our  study^  so  far  has  been  restricted  to  the 
quasi-steady"  part  of  the  mean  motion  (time  mean) . 
A tentative  shear  wave  model  for  treating  the 
nonsteady"  problem  has  been  suggested  by^  Dr. 
Ling,  and  included. 


(c)  The  flux  per  unit  gradient  is  named  the 
coefficient  of  diffusivity. 

In  connection  with  this  transport  equation,  the 
following  two  studies  are  of  particular  interest: 

(a)  Simplification  of  the  equation  to  corre- 
spond with  specific  physical  problems; 

(b)  Effect  of  turbulence  on  diffusivity,  K, 
and  hence  on  its  solution  d. 

For  lower  atmospheric  study  the  following 
simplifications  are  usually  obtained. 

1.  Vertical  stratification  and  horizontal  wind 

V=Z:  q^=iU\  so  q^-Vd  = 0 


and  equation  39  becomes 


bt  bz 


(40) 


which,  upon  putting  *S  = 0,  forms  the  basic  equa- 
tion for  usual  nonsteady,  micro-meteorology 
study; 


bt  bz 


(40a) 


2.  Steadiness  condition,  ^=0 


TURBULENT  TRANSPORT  IN  LOWER 
ATMOSPHERE 

Equation  of  Turbulent  Transport 


Equation  39  further  reduces  to 


(41) 


The  transport  of  any  physical  quantity  in  the 
atmosphere  is  governed  by  the  following  basic 
differential  equation: 

f^^^^+i■^e=V.(KV9)+S(0)  (39) 

in  which  6 denotes  the  transported  quantity,  q 
the  wind  velocity,  K the  coefficient  of  diffusivity, 
and  S the  source  function.  Physically  the  equa- 
tion states  that: 

(a)  The  total  derivative  of  quantity  6 is 
given  by  the  smn  of  divergence  of  flux 
and  rate  of  production  of  6] 

(b)  The  flux  of  6 is  in  turn  proportional  to  the 
gradient  of  0; 

673752  0—63 3 


Equation  41  actually  forms  the  basis  for  the 
first  two  sections  of  this  report  (pp.  3 and  7) 
and  for  the  present  investigation.  Upon  putting 
*^=0,  equation  41  simplifies  to 


i.e.,  the  well-known  “constant  flux”  equation  for 
steady  earth  surface  layer. 

For  the  present  report,  B will  denote  either 
carbon  dioxide  concentration  in  photosy^nthesis 
study  or  horizontal  windspeed,  u,  in  velocity  dis- 
tribution study. 

It  is  clear  from  the  above  that  our  study 
naturally  separates  into  two  regions  by  the  crop 
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height,  z=h.  For  the  outer  flow  (z^h),  S=0, 
the  governing  equations  for  both  [CO2]  and  u are 
given  by  equations  40a  and  41a;  i.e.,  the  well- 
known  atmospheric  diffusion  equation.  For 
quasi-steady  case,  we  have  simply: 


For  the  inner  flow  (z<Ch),  on  the  other  hand, 
S^^O.  In  this  case  it  has  been  shown  that  for 
quasi-steady  flow  we  have  (35) : 

[CO2I  study:  ~ (k  |^)+C(0)9+£f(2)=O;  (43) 

« study:  ^ (k  ^^+5(2)92=0.  (44) 

Turbulence  and  Dififusivity 

The  role  played  by  turbulence  in  atmospheric 
transport  takes  on  the  form  of  coefficient  of 
diffusivity,  K.  Prandtl  first  introduced  his  fa- 
mous mixing  length  theory  for  the  mechanism  of 
turbulence  transport  and  showed  that  K is 
actually  proportional  to  his  mixing  length,  1.  By 
further  postulating  that  in  the  neighborhood  of 
a solid  wall,  mixing  length,  I (and,  consequently, 
diffusivity,  K)  grows  linearly  with  distance  s from 
the  waU,  he  was  able  to  deduce  from  equation  42 
that  the  velocity  profile  for  turbulent  wall  flow 
must  follow  a logarithmic  law;  i.e.. 

This  velocity  profile  has  indeed  been  confirmed 
by  many  observations. 

In  the  last  20  years,  however,  there  has  been 
a noted  trend  among  investigators  toward  taking 
a diametrically  opposite  view  in  discarding  the 
mixing-length  hypothesis  completely  and  accept- 
ing the  logarithmic  law  as  the  basis  for  theoretical 
development,  from  which  the  “austausch  coeffi- 
cient” K can  be  shown  to  grow  linearly  with  2. 
So  far  as  the  outer  flow  is  concerned,  it  makes 
no  difference  which  point  of  view  to  take.  ’ How- 
ever, as  soon  as  we  enter  into  the  inner  flow 
region,  where  there  are  distributed  sources,  both 
approaches  face  impossible  difficulty.  On  the 
other  hand,  by  extending  PrandtPs  linear  growth 
hypothesis  instead  to  the  canopy  flow  diffusivity, 
i.e.,  K'^z,  equation  44  becomes: 

So  we  were  able  to  compute  the  theoretical  velocity 
profile  and  predict  a curvature  reversal  deep  in 


the  crop,  which  has  been  amply  verified  by  various 
field  measurements. 

This  fact  is  considered  to  be  of  great  significance 
toward  justifying  PrandtPs  mixing  length  concept 
in  general,^  and  pointing  to  the  basic  correctness 
of  our  theoretical  canopy  flow  model  in  particular, 
i.e.,  turbulent  flow  with  linearly  growing  diffusiv- 
ity, K. 

Theory  of  Average  and  Its  Applications 
1.  Theorems 

The  following  definitions  and  theorems  are  in- 
tended for  clarifying  some  basic  questions  con- 
cerning the  process  of  ensemble,  space,  and  time 
averages  employed  in  the  atmospheric  turbulence 
study. 

Definition. — The  mean  value  of  a bounded  fluctu- 
ating function  / is  defined  by: 

.mdt' . 

Theorem  1 . — It  is  evident  that 

(a)  /(Olr^o  =/(0; 

(b)  /(0|r-^o==const.; 

(c)  kt)\T  -Mt,T).  (47) 

Theorem  2. — The  mean  of  a simple  periodic 
function  retains  its  periodicity.  To,  but  the  ampli- 
tude falls  off  according  to  (T/To)~h  It  may  be 
noted  that  there  exists  a set  of  eigen-functions: 

/(^)|r=nro^O.  (48) 

Definition. — The  spectrum  of  a fluctuating  func- 
tion/ is  defined  by: 

/(„)|^=J//(0e‘“‘rf«=/(a>,  r);  (49) 

/(^:)|x=J^  j\x)e"“dx=J(k,X).  (49a) 

Theorem  3. — It  is  evident  that  /(io,oo  )=/(co), 

i.e.,  the  Fourier  Spectra. 

Definition. define  the  two  following  time 
and  space  correlation  functions: 

R{T)=^u{x,t)u{x,t-[-T)]  (50) 

2 Indeed  it  would  be  interesting  to  introduce  a parametric 
representation  of  the  “Austausch”  coefficient  as  follows: 
K'^z^  where  a is  a function  of  S.  as  8— >0.  a{S)  can 
be  determined  by  correlating  the  canopy  flow  velocity 
proflle  with  varying  crop  density. 
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t);  (51) 

where  (“)  denotes  ensemble  average. 

Theorem  4- — Ergodic  Principle;  If  a turbulent 
field  has  a stead}^  time  mean  and  a uniform  space 
mean,  then  the  field  is  ergodic ; that  is,  the  ensemble 
average  can  be  replaced  by  either  time  or  space 
average. 

In  the  present  case  this  gives 

i?(r)=Lim  = 4 f u{x,t)u{x,t-\-T)dt]  (50a) 
T->co  u'^  Jo 

i?(^)=Lim=i  f u{x,t)u{x-\-^,t)dt.  (51a) 

T->oo  u'^  Jo 

Infinite  extent  of  T ov  X cannot  be  reached, 
and  so  a real-steady,  uniform  mean  does  not 
exist.  This  state  can,  however,  be  approximated 
^ a quasi-steady  and  quasi-uniform  mean  over 
T and  X.  T denotes  the  maximum  averaging 
period  or  length  of  record  and  X the  actual  field 
extent. 

For  a turbulent  field  with  quasi-steady,  quasi- 
uniform mean,  quasi-ergodic  hypothesis  permits 
the  following  replacement: 

/(  ) =2^  r J] 

where  T is  the  finite  averaging  period,  Xn=x(tn), 
and  This  gives — 

Theorem  5. — 

= = r u{Xnt)u{Xnt  + T)dt  = Rr{T)] 

U ^ iV  1 71=1  Jo 

(53) 

■R(f)=4lwS  f u(x,,t)u(x„+(,t)dt=Ri(T). 

U ^ iV-f  71  = 1 Jo 

(54) 

Theorem  d.— For  a bounded,  fluctuating  func- 
tion f{x,  t)  with  quasi-steady,  quasi-homogeneous 
mean,  we  postulate  the  following  similarity  law: 

Br{T)^EdT);  (55) 

i.e.,  as  functions  of  averaging  period,  T,  the  space 
and  time  correlation  functions  behave  in  similar 
manner. 

The  general  behavior  of  Br  and  i?|,  as  functions 
of  T in  the  lower  atmospheric  turbulence,  is 
shown  in  reference  {24:  fig.  1.2.1,  p.  40)  and 
reference  {3:  fig.  3. 4.15,  p.  201),  respectively. 
From  the  figures  given  in  the  references,  it  is 
evident  that  both  Br  and  B^  are  monotonically 
increasing  functions  of  T,  and,  indeed. 


Lim  BrfiT)=B{T,  ^).  This  means  Bt,^{T)  is 

T-Xxi 

a faster  decreasing  function  than  Br.^  with  increas- 
ing r,  ^ 

Hence,  by  the  Birkhoff  Ergodic  Theorem  {12), 
the  applicability  of  a quasi-ergodic  hypothesis  is 
assured. 

Theorem  7. — Eddy  scale  L is  defined  by  {27)\ 

(56) 

=aj^  R,{T)dT=aL,{T)  (56a) 

from  theorem  6.  Thus,  the  recognizable  eddy 
scale  grows  with  the  interval  of  average  T. 
Theorem  8. — In  equation  56a,  we  can  put 

a=U.  (57) 

By  field  measurements,  Cramer  {4)  established  the 
following  remarkable  relationship, 

Br{T)=B^{T),^=lh,  (58) 

based  on  a 20-minute  sampling  interval  {T)  with 
U the  20-minute  mean  windspeed.  These  cor- 
relation curves  fit  approximately  the  Kolmo- 
goroff  Law  for  inertial  subrange;  i.e., 

l_i?^^2/3  (59) 

and  the  agreement  improves  with  increasing 
turbulence  level,  indicating  that  isotropy  extends 
further  into  lower  frequency  range  with  higher 
turbulence  intensity. 

2,  Application  to  atmospheric  turbulence 
analysis 

An  atmospheric  wind  record  exhibits  random 
fluctuations  containing  many  frequencies,  and 
usually  a few  dominant  ones.  The  maximum 
averaging  period,  T,  is  limited  from  above  by  the 
length  of  the  record  available. 

An  instrument  has  finite  resolution  and,  hence, 
can  only  read  a corresponding  average  with 
T>0.  The  minimum  averaging  period,  T,  is  thus 
limited  from  below  by  the  instrument  sensitivity. 
An  ideal  spectrum  V^(<^)  corresponds  to 
, T-^0.  Now  what  are  the  effects  of 

finiteness  of  T and  T on  a spectrum? 

Theorem  5 shows  that  the  maxii^um  recog- 
nizable eddy  scale  is  limited  by  T.  [Indeed 
L{T)  is  a monotonically  increasing  function  of  T.] 
Thus,  T effectively  offers  a low  frequency  (low- 
wave  number)  cutoff,  in  the  sense  that  beyond 
this  point  the  spectrum  function  defined  by 
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equation  39  only  holds  formally  and  loses  its 
physical  significance. 

Theorem  2 and  equation  49  show  that  the 
multiplying  factor  for  averaging  over  period  T 
of  a periodic  function  with  periodicity  Tq  varies 
as  (T/To)~h  Thus,  because  of  the  inevitable 
averaging  character  of  an  instrument,  the  spec- 
trum of  higher  frequency  (To<T)  or  small- 
scale  fluctuating  components  are  effectively  sup- 
pressed. In  other  words,  this  averaging  process 
removes  a large  part  of  the  time  dependence 
due  to  small-scale  eddies,  but  retains  the  time 
dependence  of  large-scale  eddies.  The  T of 
instrument  sensitivity,  therefore,  imposes  the  high 
frequency  (high-wave  number)  cutoff. 

In  relation  to  the  relevant  eddy  scale,  L{T), 
for  vertical  atmospheric  turbulent  transport  proc- 
esses, the  smaller  scale  eddies  play  the  role  of  turbu- 
lent mixing  (austausch),  while  larger  eddies  appear 
as  nonsteady  wind. 

If  there  were  no  dominant  scales  of  atmospheric 
eddies,  then  the  word  “mean”  would  have  little 
meaning.  In  the  atmospheric  surface  layer, 
however,  there  are  usually  two  natural  dominant 
scales  of  eddies : the  small-scale  ones  produced  by 
the  roughness  elements,  such  as  plant  leaf  and 
stalk;  and  the  large-scale  ones  produced  by  the 
terrain  irregularities,  like  hill  and  dale.  The 
existence  of  these  definitely  separated  ranges  of 
eddy  scale  has  been  demonstrated  by  field  record ; 
it  gives  a unique  meaning  to  the  word  “mean,” 
and,  moreover,  it  introduces  in  a natural  way  the 
concept  of  relevant  scale  and  proper  averaging 
period. 

ANALYSIS  OF  ATMOSPHERIC  WIND 
IN  THE  VICINITY  OF  PLANTS 

Background  Information  of  This  Study 

The  atmospheric  wind  consists  of  a cascade  of 
turbulent  eddies,  varying  in  size  from  a few  meters 
at  the  plant  to  macroscale  at  the  fringe  of  the 
atmosphere  {36).  The  wind  in  the  vicinity  of  a 
crop  is  influenced  by  all  scales  of  atmospheric 
disturbances.  Thus,  strictly  speaking,  no  station- 
ary time  average  exists.  However,  if  we  limit 
our  interest  to  a finite  scale  and  period,  a quasi- 
steady average  can  usually  be  found. 

The  basic  turbulent  mechanism  in  the  vicinity 
of  a crop  is  controlled  primarily  by  the  roughness 
effect  of  the  crop  alone.  In  addition,  this  mech- 
anism is  maintained  by  a driving  wind,  C/q,  that, 
within  a given  period  of  interest,  can  be  repre- 
sented by  a mean  wind  and  Fourier  spectrum  of 
oscillating  winds.  This  driving  wind  is  a function 
of  all  external  influences  other  than  that  of  the 
crop.  With  this  in  view,  it  is  possible  to  treat 
the  turbulent,  atmospheric  boundary  layer  as 
three  intermixed  states  of  flow,  namely:  (1) 
steady;  (2)  quasi-steady;  and  (3)  nonsteady. 


The  mean  driving  wind  contributes  to  the  local 
steady  state,  while  the  short-wave  spectrum  of  the 
fluctuating  wind  constitutes  its  turbulent  mech- 
anism. A sustained  local  wind  or  a long  wave  of 
1,000  meters  can  be  taken  as  quasi-steady  flow. 
A wavelength  in  the  range  of  40  to  1,000  meters  is 
considered  nonsteady.  The  oscillating  compo- 
nents in  this  range  propagate  downward.  Very 
little  is  known  about  nonsteady  turbulent-bound- 
ary-layer  flow.  Indeed,  the  types  of  nonsteady 
boundary  laj^ers  investigated  are  either  laminar- 
boundary layers  (5,  17,  21,  38)  or  turbulent- 
boundary layers  (11)  driven  by  a nonturbulent, 
oscillatory  flow  field  (potential  flow).  A bound- 
ary layer  driven  by  a turbulent  shear  flow  remains 
a subject  yet  to  be  studied. 

The  most  important  feature  that  governs  the 
turbulent  transport  in  the  vicinity  of  a crop  is  the 
steady -state-velocity  profile.  To  obtain  these  pro- 
files by  field  measurement,  the  sampling  period 
required  should  be  greater  than  the  relaxation 
period  of  the  boundary  layer.  The  relaxation 
period  will  be  defined  in  the  section  on  shear-wave 
theory.  For  flow  over  a crop,  this  relaxation  time 
is  approximately  10  seconds.  An  averaging  period 
of  30  seconds  should  smooth  out  the  high-fre- 
quency waves  due  to  turbulent  eddies.  If  the 
averaging  period  is  extended  to  cover  the  maxi- 
mum range  of  the  nonsteady  wind,  the  time  de- 
pendent factor  of  this  component  can  again  be 
eliminated.  However,  it  can  be  shown  from  field 
measurements  (fig.  10)  that  the  intensity  of  non- 
steady as  well  as  quasi-steady  winds  is,  in  general, 
beyond  the  linear  perturbation  range.  Conse- 
quently, the  result  of  a long  averaging  period  is  to 
include  nonlinear  effects  in  the  steady-state 
average. 

So  far,  the  short-period  averaging  seems  to  be 
most  logical  for  this  study,  but  this  must  be 
accompanied  by  a selective  sampling  technique 
to  eliminate  the  nonsteady  components  as  well 
as  the  nonlinear  effects  due  to  extreme  variation 
of  wind  velocities.  The  method  used  in  this 
report  embraces  the  method  of  selecting  samples 
from  periods  of  30-second  duration  preceded  and 
followed  by  approximately  the  same  time  average. 
Field  data,  obtained  by  this  means,  were  found  to 
be  quite  consistent.  This,  in  a way,  verified  the 
basic  concept  of  selective  short-time  average. 

Wind-profile  measurements  over  a wheatfield 
and  a cornfield  were  conducted  with  multiple 
small-cup  anemometers  of  9-cm.  diameter,  in 
conjunction  with  Hasting  thermocouple-anemom- 
eters.^ The  cup  anemometers  were  appropri- 
ately mounted  at  various  heights  within  and 
above  the  crop.  Thermocouple  anemometers 
were  used  within  the  crop  where  the  windspeed 
was  considered  too  low  for  cup  anemometers. 

* The  use  of  this  or  other  patented  equipment  in  this 
study  does  not  imply  approval  of  the  product  to  the 
exclusion  of  others  that  may  also  be  suitable. 
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Figuee  10. — Simultaneous  horizontal-wind  recording  at  two  different  elevations  above  a corn  crop. 


Simultaneous,  integrated  readings  from  these 
anemometers  were  recorded  at  30-second  interTals. 
The  average  wind  profiles  were  obtained  from 
these  measurements.  The  selected  samples  rep- 
resenting quasi-steady  wind  profiles  for  a wheat 
field  and  a cornfield  are  shown  in  figures  11  and  12, 
respectively.  The  range  of  wind  profiles  pre- 
sented here  is  limited  to  the  most  common  pre- 
vailing winds  encountered  in  a typical  field. 
Data  for  extremely  low  and  high  winds  are  in- 
complete, and,  therefore,  not  presented. 

The  wheatfield  was  situated  on  top  of  Mount 
Pleasant  unobstructed  on  all  sides,  whereas  the 
cornfield  was  situated  in  Ellis  Hollow.  (See 
last  section,  p.  37.)  The  wind  profile  character- 
istics, by  virtue  of  the  sampling  technique  em- 
plo3’-ed,  are  expected  to  be  independent  of  the 
geographical  terrain. 

Most  of  the  field  data  and  the  instrumentation 
work  used  in  this  report  were  carried  out  with  the 
assistance  of  J.  Stoller,  Xew  York  State  College 


of  Agriculture,  under  the  direction  of  E.  R.  Lemon, 
ARS,  USDA.  A considerable  amount  of  field 
data  had  been  compiled  during  the  summer 
months  of  1960,  but  only  a small  part  of  the 
information  was  analyzed  and  presented  herein,  as 
a consequence  of  limited  time  and  effort. 

The  general  theory'  and  experimental  results 
of  steady  state  and  quasi-steady  state  winds  in  the 
vicinity  of  a crop  are  presented  in  the  next 
subsection. 

Boundary-Layer  Profiles 

1.  Logarithmic  wind- profile  over  a crop 

The  logarithmic-boundary-layer  theory  of 
Prandtl  {23)  applies  very  well  to  the  steady- 
bound  arj^-l  ay  er  flow  above  a crop.  In  the  original 
theoretical  model,  it  was  assumed  that  the  mixing 
length  or  the  equivalent  eddy  viscosity  of  turbulent 
mixing  grows  linearly  with  the  height.  This 
assumption  is  now  extended  all  the  way  from 
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Figure  11. — Wind  profiles  over  a wheat  crop,  1960. 


Z 

(meters) 


Figure  12. — Wind  profiles  over  a corn  crop,  1960. 

and  incorporating  equations  61  and  62  in  equation 
60,  we  obtain  the  well-known  logarithmic- 
boundary-layer equation  through  integration  of 
equation  60. 


within  the  crop  through  the  inner  boundary  layer 
up  into  the  outer  turbulent-flow  layer  (21).  By 
virtue  of  our  stratification  hypothesis;  i.e.,  van- 
ishing pressure  gradient  we  have 


i.e.,  the  flow  is  maintained  through  the  shearing 
action  of  eddy  viscosity.  Now  this  leads  to 

r = e(Z)  ^=constant.  (60) 


U(Z,Uo)_l.  Z-d(Uo) 
U,iUo)  ZoiUo)  ’ 


zyh 


(63) 


where 


^=0.42  Von  Karman’s  constant; 

Z— height  above  ground  level; 

C7o=Diean  wind  velocity  at  10-meter  reference 
level ; 

^=peak  plant  height,  defined  as  1.1  times 
the  mean  plant  height,  h; 
d = zero-plane  displacement,  a function  of  Uo ; 
Zo= roughness  length,  a function  of  Uo; 
p=mass  density  of  air. 


Introducing  the  friction  velocity  as 

= (61) 

the  eddy  viscosity  €(Z),  as 

e{Z)  = pU^k{Z—d),  (62) 


The  constants,  U*,  Zo,  and  d are  called  the  char- 
acteristic parameters  of  a steady-state-wind  profile. 

In  order  to  provide  a convenient  reference  level 
in  place  of  a boundary  thickness,  the  conventional 
height  of  10  meters,  Zio,  was  chosen  as  the  nominal 
limit  of  the  surface  layer.  Ratios  formed  with 
Zjo  are  to  be  taken  as  of  approximate  order  of 
magnitudes. 
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The  three  wind-profile  parameters,  d,  and  Zo, 
can  be  determined  from  the  measured  \vind  pro- 
files, as  shown  in  figures  11  and  12,  and  they  are 
obtained  through  the  following  equations: 


and 


U^=k 


U,-U2 


In 


iZ^-d) 

{Z2-d) 


z„= 


Z,-d 

^kUijU*  ■ 


(65) 


(66) 


These  equations  are  derived  directly  from  the 
basic  logarithmic  equation  63.  The  values  of  these 
parameters  as  a function  of  reference  mnd,  Uq,  for 
a wheatfield  and  a cornfield  are  plotted  in  figures 
13  and  14,  respectively.  It  is  to  be  noted  that 
values  given  in  the  above  figures  are  for  a typical 
wheatfield  and  a cornfield  having  an  approximate 
plant  density  of  60  plants  and  0.5  plant  per  square 
foot  of  ground,  respectively. 

The  value  of  kinematic-eddy  viscosity.  Km, 
known  also  as  the  coefficient  of  eddy  diffusivity, 
is  of  basic  interest  in  the  study  of  diffusion  of 
CO2,  water  vapor,  and  heat.  Direct  measure- 
ment of  this  coefficient  was  planned,  but  only 
crude  measurements  were  possible  at  the  time 
of  this  report.  In  general,  they  check  well  with 
those  obtained  indirectly  from  the  wind-profile 
parameters.  The  equation  for  eddy  diffusivity 
is,  from  equation  62, 


K„=^=U^{Uo)k[Z-d(Uo)].  (67) 

P 


Values  of  and  d can  be  obtained  from  figures 
13  and  14  for  a given  Uo  and  field  condition. 

The  drag  coefficient  defined  as 

2 

is  also  plotted  in  figures  13  and  14. 

Comparing  the  behavior  of  flow  parameters 
between  a wheat  and  a cornfield,  it  is  noticed 
that  the  zero-plane  displacement,  d,  decreases 
with  an  increase  of  Uq;  this  characteristic  is 
controlled  primarily  by  the  mechanism  of  the 
canopy  flow.  The  values  of  friction  velocity, 
are  of  about  the  same  magnitude  and  char- 
acteristic for  both  the  wheatfield  and  the  cornfield. 
This  indicates  that  the  characteristics  of  turbu- 
lence above  either  crop  are  of  approximately  the 
same  scale  and  intensity.  The  value  of  roughness 


parameter,  Zo,  behaves  differently  for  the  two 
crops.  For  the  wheatfield,  Zo  tends  to  decrease 
with  Uo,  whereas  the  reverse  is  true  for  the 
cornfield.  It  is  postulated  that  the  behavior  of 
Zo  depends  on  the  fiexibility  of  the  plant  as  well 
as  plant  density.  The  wheat  plant  is  more 
flexible  and  tends  to  bend  and  “streamline” 
itself  under  a high  wind,  but  the  corn  is  relatively 
stiff  and  does  not  “streamline”  itself  as  readily. 
The  difference  in  effective  plant  density,  or  leaf 
area  density  between  wheat  and  corn,  exercises 
a strong  influence  on  the  canopy  flow  that  in- 
directly affects  the  characteristic  of  Zo.  The 
different  dependence  of  Cd  on  Uq  between  a 
wheatfield  and  a cornfield  is  similar  to  the  varia- 
tion of  Zp  Avith  Uq.  Therefore,  the  same  tentative 
explanation  regarding  plant  flexibility  is  applicable. 

2.  Canopy  How  within  a crop 

The  basic  analytical  model  for  a canopy  flow 
has  been  mentioned  in  the  introductory  section 
of  this  report.  Basically,  it  is  assumed  that  the 
flow  within  the  crop  is  maintained  by  the  action 
of  eddy  viscosity,  and  that  the  mixing  length 
theory,  which  is  applied  to  the  outer  flow,  is 
also  valid  in  this  zone.  The  differential  equation 
governing  the  flow  is, 

(69) 

where 

^4(Z)=total  vertical  plane  area  of  plants/unit 
volume,  a function  of  Z;  and 

(71,= drag  coefficient  of  plant. 

The  term  on  the  right  represents  the  drag  force 
offered  by  leaves  and  stalks  of  the  crop.  Within 
the  flow  range  considered,  C'd  can  be  taken  as 
independent  of  Keynolds  number,  and  the  drag 
force  is  proportional  to  the  square  of  the  velocit3\ 
Introducing  equation  60  into  equation  69, 
we  have 

Let 
— Z 

Z=^»  normalized  height  variable,  and 
fi 

U=^’  normalized  wind,  where  /7ft=wind  at 
peak-plant  height, 

and  introducing  these  new  variables  in  equation 
70;  we  have 


-=+«(2) 

J dZ  dZ  ^ 

Z<1.  (71) 
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Figure  13. — Parameters  for  steady-state  wind  profiles  for  winter  wheatfield. 
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Figl're  14. — Parameters  for  steady-state  wind  profiles  for  cornfield. 
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In  order  to  define  statistically  the  value  of 
peak-plant  height,  some  sampling  data  on  the 
probability  density  distribution  of  corn  height, 
P(A),  were  taken  and  plotted  in  figure  15.  The 
peak-plant  height  has  a value  of  about  1.2  times 
the  mean  height  of  plants.  However,  instead 
of  taking  this  figure,  a value  of  1 . 1 times  the  mean 
height  of  plants  is  taken  as  the  definition  of  h. 
It  was  assumed  that  this  definition  is  also  appli- 
cable to  wheat. 

The  area  distribution  function,  A{Z),  can  be 
obtained  statistically  for  a given  crop.  However, 
for  simplicity,  it  was  assumed  a priori  that  the 
average  vertical  projection  area  of  a single  plant 
(with  the  vertical  plane  making  an  angle  of  45° 
between  the  planes  of  maximum  and  minimum 
projected  area)  has  an  area  distribution  that  is 
uniform  from  the  root  of  a plant  to  85  percent  of 
its  height;  from  that  point  on,  the  area  reduces 
uniformly  with  height  to  zero  at  the  tip.  By 
combining  this  standard  area  function  with  the 
plant  height  distribution  function,  the  statistical 


area  distribution  function,  A{Z),  is  obtained  and 
plotted  as  a ratio  of  area  density,  Aq,  at  the 
midplant  height.  This  function  is  shown  with 
“x”  marks  in  figure  15.  The  area  distribution 
function  again,  for  the  sake  of  simplicity,  was 
approximated  by  a straight  line  shown  "in  the 
above  figure.  This  simplified  area  function  gives 
zero  area  at  1.1  times  the  mean-plant  height, 
and  hence  this  height  was  used  for  the  definition 
of  peak-plant  height. 

If  we  introduce  the  linear  growth  of  eddy  vis- 
cosity with  height, 

where  is  the  eddy  viscosity  at  height,  h,  and  with 
the  above  simplified  function  of  A{Z),  equation  71 
is  written  for  two  zones : 

for  0<Z<0.73 
A{Z)=Aq,  constant 


P (h)% 

0 20  40  60  80  100 
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Figure  15. — Plant  area  and  height  distribution  functions 
for  corn. 


(72) 


for  0.73<Z<1.Q 


{Z-o.7Z)-\u\ 


(73) 


where  S'  is  the  dimensionless  parameter  of  the 
differential  equation.  It  is  to  be  called  the  shape 
factor  of  canopy  flow,  namely, 

_phHJ  }iAqCd 

Equations  72  and  73  are  nonlinear,  differential 
equations  of  the  second  order.  They  are,  how- 
ever, unique  and  can  be  solved  by  numerical 
integration  for  the  following  boundary  conditions : 

7/z=i  = l ; 

Uz=o=0. 

Relaxation  method  (28)  was  used  to  obtain  the 
characteristic  canopy-wind  profiles  for  various 
values  of  shape  factors.  The  ^result,  plotted  in 
figure  16,  covers  the  range  of  S from  5 to  1,000. 

In  all  these  cases,  the  ground  shear  w^ 
neglected  in  comparison  to  the  leaf  drag.  ]^r  S 
<5,  ground  effect  comes  into  play.  When  S-^0, 
ground  shear  alone  controls  the  flow,  and  we 
obtain  the  normal  logarithmic  wind  profile. 
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Figure  16. — Canopy  wind  profiles  for  various  shape 
factors  S. 


3.  Combination  of  canopy  and  logarithmic- 
wind  profile 

The  canopy  flow  and  the  logarithmic-profile  flow 
above  are  actually  one  continuous  system.  The 
matching  point  between  the  two  flow  fields  is  at 
the  peak-plant  height,  h,  where  the  wind  velocity 
and  the  velocity  gradient  must  match.  This 
means  that,  for  a given  logarithmic-wind  profile, 
there  will  be  a canopy  profile  of  a given  shape 
factor  that  will  match  it  at  A.  The  value  of  en  in 
equation  74  must  also  match  the  value  of  e(Z);, 
determined  by  equation  62  of  the  logarithmic-wind 
profile. 

The  velocity  gradient  ^ for  the  canopy 
dZ  ^=1 

flow  is  plotted  in  figure  17  as  a function  of  shape 
factor  S.  From  that  plot,  the  matching  shap^ 
factors  for  the  wind  profiles,  shown  in  figures  11 
and  12,  are  plotted  in  figures  13  and  14,  respec- 
tively. The  trend  of  shape  factor  for  the  cornfield 
and  the  wheatfield  turns  out  to  be  the  same;  that 
is,  it  tends  to  decrease  with  an  increase  of  reference 
wind,  Uq.  The  corresponding  values  of  wheat 
are  about  10  times  those  of  corn.  The  reason 
for  this  difference  is  that  the  effective  area  density, 
Aq,  of  the  wheatfield  is  much  higher  than  that  of 
the  cornfield. 

In  figure  16  the  measured-canopy  flow  profiles 
for  corn  and  for  wheat  are  plotted.  They  fit  well 
with  the  computed  wind  profiles.  This  validates 


iQ| 


Figure  17. — Velocity  gardients  at  peak-plant  height  as 
a function  of  shape  factor  S. 

the  assumption  of  the  turbulent  mixing-length 
concept  used  in  the  analytical  model. 

For  both  the  wheatfield  and  the  cornfield,  the 
ratio  of  eddy  viscosity  at  peak-plant  height  to 
the  viscosity  of  air  varies  from  1,000  to  8,000 
within  the  range  studied.  The  assumption  of  a 
fully  turbulent  model,  and,  consequently,  the 
application  of  mixing-length  theory  to  the  study 
of  air  motion  in  the  vicinity  of  a crop,  are,  there- 
fore, valid. 

Shear-Wave  Theory 

During  the  course  of  this  study,  a very  crude 
shear-wave  theory  was  proposed  to  study  one 
phase  of  the  nonsteady  state  phenomena.  This 
idea  was  generated  by  field  measurements  of  the 
wind  profiles  when  it  was  observed  that  wind 
components  at  different  heights  exhibited  phase 
lags,  thereby  suggesting  a propagating  shear- wave. 
Figure  10  shows  two  simultaneous  recordings  of 
the  horizontal  wind  at  two  vertical  positions 
Z=560  cm.  and  Z=300  cm.— within  a cornfield. 
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The  lines  containing  high  frequencies  are  instan- 
taneous recordings  from  anemometers  with  a 
response  time  of  less  than  one  second.  The 
second  set  of  lines  in  long  dashes  represents  the 
running  average  of  15  seconds,  and  the  third  set 
represents  a running  average  of  60  seconds.  For 
the  60-second  running  average,  all  wavelengths 
under  60  seconds  are  averaged  out,  and  only  the 
wave  spectrum  greater  than  60  seconds  remains. 
These  long-period  waves  vary  nearly  in  phase 
with  each  other,  and  are  defined  as  quasi-steady 
flow.  The  15-second  average  line,  which  oscillates 
about  the  60-second  average  line,  represents  the 
spectral  components  between  15-second  and  60- 
second  waves. 

These  waves  are  seldom  in  phase  between  the 
two  levels  under  observation.  The  top  wave 
may  either  lead  or  lag  behind  the  lower  waves. 
It  is  not  easy  to  differentiate  which  is  leading 
or  lagging  when  waves  of  multiple  periods  are 
intermixed.  However,  when  there  is  a sudden 
increase  in  wind  intensity,  the  wind  on  the  upper 
level  generally  can  be  observed  to  be  leading 
that  of  the  lower  level.  This  type  of  oscillating 
wind  is  identified  as  nonsteady  wind. 

The  differential  equation  for  shear-waves  prop- 
agating in  the  Z direction  can  be  obtained 
through  the  modified  Navier-Stokes  equation  by 
replacing  eddy  viscosity  for  molecular  viscosity 
and  neglecting  pressure  gradients  and  velocity 
gradients  in  the  x and  y directions. 


P 


r (7\ 

dt  ’ bZ\  dZ 


(75) 


Differentiating  equation  75  with  respect  to  Z 
and  multiplying  each  side  by  €(Z),  we  have 

a w €(Z)  av 
dZ  p dZ^ 
or 

%=K(Z)  (76) 


The  eddy  diffusivity  i^(Z)  may  be  taken  as 
linearly  proportional  to  Z;  however,  this  leads  to 
a rather  complex  solution  for  the  shear  wave.  To 
simplify  the  situation,  it  is  expedient  to  consider 
that  K{Z)  is  a constant,  with  an  average  value 
of  Kq  within  the  range  considered.  Equation  76 
is  then  reduced  to 


av 

az^' 


(77) 


with  the  boundary  conditions: 


Z=0 


t=tq-\-t'  sin  (cor). 

The  steady  state  solution  for  equation  77  is 
similar  to  the  well-known,  periodic-heat-transfer 
problem  of  a fiat  slab  (2). 


where 


r(Z,  t)  = To+ 


t'  cosh  oi'Z{l-[-i) 
cosh  co'Zio(l+^)  ^ 


(78) 


When  the  frequency  factor  co'Zjo  is  less  than 
0.5,  the  absolute  magnitude  of  oscillating  shear, 
t'(Z),  differs  only  2 percent  at  Z=0  from  r'  at 
Zio.  The  corresponding  phase  lag  at  these  levels 
is  14°.  That  is,  when  the  frequency  factor  is 
less  than  0.5,  constant-shear  condition  can  be 
assumed. 

It  can  be  shown  that  the  above  shear  wave  will 
propagate  downward  with  a speed  of 


Oi 

dr  cosh  o}'Z{l-\-i)  "I  (79) 
dZ  cosh  co'Zio(l-f  f)  J 

The  speed  of  wave  propagation  is  not  constant 
with  respect  to  height.  However,  it  can  be 
shown  that  for  a good  approximation,  Vr  may  be 
taken  as 

^>0.2.  (80) 
■^10 

The  shear-wave  length  in  the  Z direction  can 
be  defined  as 


X 

=2Tr-\ 

(81) 

^ CO 

\ 

/ CO 

The ‘ratio  of  Zio  to  X 

r is 

— — — Z 

/ “ - 

1 

o 

(82) 

27r 

For  the  limiting 

case 

of  constant 

shear, 

dr/dZ—O,  where  Zjoco'  is  less  than  0.5,  the  ratio  of 

^"<^=0.08. 

At  aTT 

That  is,  when  the  shear-wave  length  is  more  than 
12  times  the  reference  height,  Zio,  constant-shear 
condition  exists  (quasi-steady  state) . The  shear- 
wave  length  required  for  this  state  of  flow  is 

Xr>  7^=125  meters. 

U.Uo 


and  r at  reference  height  Zjo  oscillates  about  a 
mean  shear  tq  as 


From  equation  81  the  frequency  of  oscillation,  /, 
for  quasi-steady  condition  can  be  obtained. 
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(83) 


For  a typical  condition  where  Uo=500  cm. /sec., 
X^  = 12,500  cm.,  and  i^o=50,000  cm.^sec.,  the 
corresponding  value  of / from  equation  83  is 


/= 


47r50,000 

12,5002 


=4X10- 


cycles/sec. 


Under  the  conditions  when  the  time  and  space 
coordinates  are  interchangeable,  the  equivalent 
velocity  wavelength  can  be  defined  as 


500 

4X10-3 


= 125,000  cm.  = 1,250  meters. 


Hence,  from  the  very  crude  shear -wave  theorj^,  an 
approximate  boundary  value  for  quasi-steady 
state  condition  is  defined.  That  is,  an}"  wave- 
length greater  than  1,000  meters  may  be  considered 
as  a quasi-steady  wind. 

The  wind  spectrum  can  now  be  divided  into 
the  following  classification  of  edd}"  scales: 


Scale  Ratio  Eddy  Classification 

^<4 Turbulent  eddies. 

^10 


4<^<100 Nonsteady  wind. 

•^10 


100<^ Quasi-steady  wind. 

■^10 

The  relaxation  time,  U,  of  a boundary  layer, 
referred  to  at  the  beginning  of  this  section,  is 
defined  as  the  time  for  the  shear  wave  to  propagate 
from  the  reference  height  down  to  the  plant. 

(84) 

For  the  typical  example  as  given  above,  the  value 
for  tr  is  14  seconds.  Winds  of  higher  frequenc}"  or 
shorter  wavelength  will,  of  course,  give  shorter 
periods.  The  phase-time  relationship  of  waves 
of  30  seconds  duration,  as  shown  in  figure  10  (in 
long  dashes),,  is  about  4 seconds  between  the  two 
recording  levels.  Equation  84  gives  a relaxation 
time  of  about  the  same  magnitude.  This  checks 
well  with  the  experimental  data. 


ON  THE  DIFFUSION  OF  HEAT  AND 
MASS  IN  A TURBULENT  SHEAR 
FLOW 


By  D.  A.  Spence 

The  feature  that  distinguishes  turbulence  from 
laminar  flow  is  its  randomness:  In  any  given 


case  it  is  possible  to  give,  at  the  very  best,  only  a 
statistical  description  of  the  way  in  which  the 
distributions  of  velocity,  pressure,  temperature, 
and  density  vary  about  their  means.  Fortu- 
nately, such  a description  is  usually  good  enough, 
since  one  is  only  interested  in  the  behavior  of 
properties  averaged  over  a time  that  is  long 
compared  with  that  of  a turbulent  fluctuation. 
Exactly  as  in  statistical  mechanics,  one  is  only 
interested  in  the  mean  properties  of  an  ensemble, 
and  one  not  only  could  not  but  would  not  want 
to  say  how  an  individual  particle  behaves. 

The  analogy  with  statistical  mechanics  is  a 
natural  one,  and  led  Prandtl  to  introduce  the 
notion  of  a mixing  length  as  a counterpart  of  the 
mean  free  path,  but  it  is  nowadays  generally 
agreed  that  although  this  concept  has  dimensional 
validity,  it  is  a great  oversimplification  of  the 
transfer  processes  at  work  in  a turbulent  flow. 
For  this  reason,  one  is  on  dangerous  ground  in 
trying  to  predict  the  results  of  a new  experi- 
ment— for  example,  in  appljfing  suction  or  blo^ving 
to  a turbulent  boundary  layer — by  means  of 
mixing  length  theory,  even  though  it  can  be  used 
very  successfully  to  rationalize  some  well-known 
observations,  the  most  notable  being  the  loga- 
rithmic law  of  the  wall. 

Turbulence  presents  a much  more  difiicult 
problem  than  the  mean  motion  of  molecules  in 
a gas,  however,  since  the  motion  of  one  fluid 
particle  affects  that  of  every  other,  and  the 
simplification  of  “molecular  chaos”  is  not  avail- 
able. Moreover,  the  governing  equations  are 
nonlinear,  and  the  fluctuating  motion  is  always 
three-dimensional,  even  in  a one — or  two — di- 
mensional mean  flow.  For  this  reason,  substan- 
tial progress  in  predicting  the  development  of  a 
field  of  turbulence  from  a statistically  known 
initial  state  has  only  been  possible  in  the  simplest 
case : that  of  decaying  homogeneous  isotropic 
turbulence.  Definitive  accounts  of  this  subject 
are  given  by  Batchelor  {1)  and  Lin  {19),  and 
mention  should  also  be  made  of  the  review  by 
Liepmann  {16). 

The  primary  aim  of  workers  in  this  field  is  to 
predict  the  form  and  time  variation  of  the  turbu- 
lent energy  spectrum.  In  this  way  they  obtain 
some  insight  into  the  mechanisms  by  which  energy 
is  transferred  by  the  inertia  terms  in  the  equations 
of  motion  from  large  eddies  with  low-wave 
numbers  through  an  intermediate  range  to  the 
high-wave  numbers  at  which  dissipation  takes 
place  by  viscous  action.  The  well-confirmed 
existence  of  an  inertial  subrange  in  which  the 
spectrum  depends  only  on  the  rate  of  dissipation 
and  the  wave  number,  k — and  is,  therefore,  found 
by  a simple  dimensional  argument  to  vary  as 
^-5/3 — shows  that  the  process  is  a continuous  one; 
i.e.,  that  energy  is  passed  successively  down  the 
whole  wave-number  scale  without  bypassing 
any  part  of  the  range.  To  make  further  prog- 
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ress,  however,  some  physical  hypothesis — for 
example,  that  the  velocity  correlations  of  a certain 
order  are  those  given  by  joint  normal  distribu- 
tions— must  be  made,  and  nothing  completely 
satisfactory  seems  to  have  emerged,  so  the  theory 
of  homogeneous  turbulence  is,  at  present,  some- 
thing of  a standstill. 

In  the  meantime,  however,  the  same  methods 
have  been  used  to  acquire  a considerable  insight 
into  the  qualitative  nature  of  turbulence  in  the 
important  class  of  flows  in  which  the  mean  motion 
is  sheared,  as  between  fixed  boundaries  in  a pipe 
or  channel,  free  boundaries,  or  in  wakes  or  jets, 
and  one  fixed  and  one  free  in  a boundary  layer. 
A comprehensive  survey  of  these  flows  was  given 
by  Townsend  (37)  in  1956.  From  his  work  a 
fairly  definite  picture  has  emerged  of  the  way 
eddies  compare  in  size  with  the  width  of  the  layer 
extract  energy  from  the  mean  flow  and  hand  it 
down  through  eddies  of  progressively  smaller  size 
to  those  in  which  it  is  ultimately  dissipated.  (The 
latter  are  not,  as  was  once  thought,  necessarily 
isotropic  even  in  a shear  flow.)  The  net  rate  of 
production  of  turbulent  energy  per  unit  mass  is 
8t 

u —)  the  product  of  mean  velocity  and  shear  stress 

gradient,  and  if  it  were  possible  by  considering 
the  structure  of  the  turbulence  to  say  how  this 
quantity  depends  on  the  local  velocity  and  pres- 
sure gradients,  the  calculation  of  mean  flow  as  a 
shear  layer  would  be  fully  determinate;  a further 
relation  being  provided  by  the  Navier-Stokes  equa- 
tions. But  Townsend  and  his  coworkers  are  still 
quite  a long  way  from  being  able  to  make  any 
quantitative  predictions  of  the  transport  properties 
in  terms  of  local  mean  quantities.  (“Transport” 
is  used  so  as  to  include  heat  and  mass  diffusion 
into  shear  force;  i.e.,  momentum  transport,  in  the 
discussion.) 

Actually,  it  is  not  at  all  clear  that  a representa- 
tion in  terms  of  local  quantities  is  physically 
meaningful,  even  though  it  is  always  formally  pos- 
sible by  means  of  an  eddy  diffusion  coefficient. 
Transport  of  any  quantity  is  brought  about  to 
different  extents  by  eddies  of  different  sizes.  In- 
tuitively one  has  a picture  of  a lump  of  fluid 
carrying  with  it  some  property  having  the  mean 
value  appropriate  to  its  initial  surroundings  and 
being  convected  by  an  eddy  to  a point  where  the 
mean  properties  are  different,  and  where  it  then 
shares  the  transported  property  with  the  surround- 
ing fluid.  The  distance  traveled  by  the  lump  is 
then  the  mixing  length,  and  the  flux  of  the  trans- 
portable property  should  be  proportional  on  this 
arg-ument  to  the  mean  gradient  over  the  length. 

But  lumps  converted  by  different  eddies  will 
travel  different  distances.  At  least  in  the  case  of 
the  layer  eddies,  an  exchange  coefficient  based  on 
the  local  gradient  at  one  point  in  the  shear  layer 
does  not  seem  appropriate,  both  because  the  eddies 
are  comparable  in  width  to  the  shear  layer  and 


because  the  exchange  process  takes  place  con- 
tinuously by  molecular  action  wherever  there  is 
a local  concentration  gradient,  not  waiting  for  the 
lump  to  reach  its  final  position.  In  fact,  Town- 
send’s experiments  (see  Hinze  {9),  pp.  288-289) 
indicate  that  momentum  transport  is  primarily 
due  to  gradient-type  diffusion  by  small-scale  tur- 
bulence (presumably  because  the  u and  v compo- 
nents in  the  large  eddies  tend  to  be  uncorrelated) , 
whereas  heat  is  transported  both  by  gradient  dif- 
fusion and  by  the  Wlk  mechanisms  associated 
with  large  eddies. 

Thus,  as  pointed  out  by  Ligh thill  {18)  in  review- 
ing Townsend’s  book,  it  is  physically  quite  signifi- 
cant to  define  an  eddy  viscosity,  e,  by 


— pu'v'  = e 


w 


and  to  expect  e to  depend  on  local  quantities.  To 
this  is  added  the  viscous  contribution  to  shear 
stress,  giving  altogether 


T=(/x-f  e) 


du 


and  the  analogy  between  the  gradient  diffusion  due 
to  eddy  and  molecular  motions  is  clear.  By 
analogy  with  the  first  equation,  an  eddy  con- 
ductivity K and  an  eddy  diffusion  coefficient  D 
(for  a primary  mixture)  may  be  defined  in  terms 
of  the  fluxes  of  enthalpy,  h,  and  mass,  denoted  by 
concentration,  c,  of  one  species,  due  to  the  turbu- 
lent motion  by 


—ph'v'=K 


by 


—pc'v'=D 


by 


The  corresponding  molecular  transport  coeffi- 
cients are  denoted  by  k and  D,  say.  (Note  that 
the  specific  heat  Cp  has  been  absorbed  in  the 
definitions  of  k and  K.)  The  kinetic  theory  for  a 
dilute  gas  at  room  temperatures  shows  that  the 
ratios 

Pr=^,  Le=^,  Sc=-^ 

k k pD 

are  all  of  order  unity  (between  0.7  and  1.4,  roughly) 
and  one  is,  therefore,  able  to  treat  the  diffusion  of 
heat  or  mass  by  analogy  with  that  of  momentum. 
For  example,  Reynold’s  analogy,  which  holds 
exactly  for  Pr=\,  states  that,  when  suitably  non- 
dimensionalized,  skin  friction  is  equal  to  heat 
transfer  for  a boundary  layer. 

For  turbulent  flows  one  naturally  seeks  to  use 
the  same  type  of  analogy,  even  though  for  the 
reasons  discussed,  it  rests  on  less  certain  physical 
grounds.  Nevertheless,  it  is  possible  to  give  quite 
good  quantitative  answers  to  practical  problems 
in  chemical  engineering,  micro-meteorology,  etc.. 
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by  relating  the  fluxes  of  heat  and  matter  to 
the  shear  stress  by  means  of  the  corresponding 
parameters 


{Pr) 


m 


which,  in  view  of  the  broadly  analogous  exchange 
mechanisms,  shoifld  all  be  somewhere  near  unity. 
General  surveys  of  this  approach  have  been  given 
by  Lees  {13)  and  by  Spalding  {30),  the  latter  being 
particularly  concerned  with  combustion  problems. 
The  usual  practice  is  to  give  the  ratios  constant 
values  in  any  given  shear  flow.  In  this  way  the 
present  author  {31)  was  able  to  relate  the  distribu- 
tion of  enthalpy  to  that  of  shear  stress,  both  being 
regarded  as  functions  of  velocity,  in  terms  of  an 
arbitrary  value  a for  the  turbulent  Prandtl  num- 
ber. If  the  energy  equation  was  integrated,  it  was 
then  possible  to  express  both  the  recovery  and 
Reynold’s  analogy-factors  in  terms  of  a,  and  by 
taking  a=0.85,  which  is  midway  between  unity 
and  (Pr)  moucuiar,  good  agreement  with  the  known 
experimental  results  for  air  was  found. 

In  view  of  the  differences  in  transport  mech- 
anisms for  heat  and  momentum,  it  would  be 
surprising  if  Pr turb  actually  were  constant  across 
the  layer  and  the  above  value  should  probably 
be  looked  on  as  a weighted  means.  The  transport 
must  depend  on  the  eddy  structure,  and  this 
varies  across  the  layer.  For  the  region  where 
small  eddies  predominate  (close  to  the  sublayer), 
gradient  diffusion  would  probably  be  the  main 
mechanism,  and  one  might  expect  to  And  (Pr)turb 
~(P^) molecular,  but  near  the  outer  edge  of  a shear 
laj^er  large  eddies  are  predominant  and  transport 
is  primarily  by  bulk  motions,  so  it  seems  more 
probable  that  the  Prandtl  number  should  then 
approach  unity.  The  integration  of  reference 
{31)  can  actual!}^  be  carried  out  formally  for  an 
arbitrary  variation  of  (Pr)turb,  as  is  shown  in  the 
notation  below,  giving  a result  of  some  generality 
for  zero-pressure  gradient  flows.  Like  most  such 
results,  they  are  not  suitable  for  application 
without  some  physical  hypothesis.  In  this 
case,  the  dependence  of  (Pr)  turb  on  the  shear  stress 
may  be  a suitable  representation  of  the  local  edd}’ 
structure. 

The  same  methods  can  be  used  to  estimate  the 
diffusion  of  mass  across  a shear  flow  in  terms  of 
shear  stress  and  velocity  in  terms  of  an  arbitrary 
Schmidt  number,  and,  in  the  absence  of  definite 
experimental  information,  one  would  take  a value 
somewhere  between  unity  and  the  molecular 
value  as  the  most  likely.  Meteorologists  usually 
use  mixing  lengths  having  some  prescribed- 
height  dependence  to  obtain  the  distribution  of 
velocity  and  shear  stress,  as  detailed  by  Priestley 
{25).  When  this  is  known,  one  can  proceed 
exactly  as  in  the  heat  diffusion  case  to  calculate 


the  flux  of  a foreign  gas  through  the  atmosphere 
due  to  turbulent  mixing. 

Integration  oj  the  energy  equation  in  a zero- 
pressure  gradient,  with  variable  turbulent  Prandtl 
number. — The  energy  equation  with  the  usual 
notation  is 


I / I ^h  dq  . du 


and  if  h = h{u),  this  can  be  written. 


1 ^ ^_1  dq  . 

T du  du  T du' 

The  heat  flux, 

I 1 

<1=(><:+K)^y=hn^< 
say,  and  likewise 

du 


If  we  write  ^=a,  then  a is  a function  of  the 
kes 


P^turb— 


Meg— M 
^eg  k 


In  terms  of  a the  energy  equation  is 


(The  wall  stress  has  been  brought  in  for  con- 
venience.) 

Denote  by  the  Stieltjes  integral 


Then,  if  we  set  a=a  at  the  waU  u=0,  the 
result  of  integrating  the  energy  equation  twice  is 

h-K=^-{f)  rj^T\-*du 

(T\dU/y,jQ  \t,c/ 

do  \r w/  do  \r w/ 

If  Q!= constant,  this  reduces  to  the  usual  Crocco 
integral. 

It  can  also  be  integrated  in  a fairly  straight- 
forward manner  if  a is  assumed  to  be  a simple 

function  of  — • For  instance,  it  might  be  physically 

r w . 

realistic  to  assume  a linear  dependence,  by  setting 

q:=:1  — (1  — 0-) 

which  goes  to  unity  at  the  free  boundary  and  to 
cr  at  the  wall. 


A METHOD  FOR  THE  COMPUTATION  OF  LOGARITHMIC  WIND 
PROFILE  PARAMETERS  AND  THEIR  STANDARD  ERRORS 


By  WiNTON  Covey 


The  logarithmic  wind  profile  equation  may  be 
written: 

y=F,  In  (^)-  (85) 

Here,  2 is  height  above  the  ground,  Fis  mean  wind- 
speed  over  a short  (e.g.,  15-minute)  time-interval. 
Equation  85  is  observed  to  fit  observations  well 
when: 

(1)  The  anemometry  is  good, 

(2)  Thermal  stratification  is  adiabatic, 

(3)  The  measuring  site  is  amid  a uniform 

patch  of  ground, 

(4)  The  lowest  anemometer  is  not  too  low, 

and 

(5)  The  highest  anemometer  is  not  too  high. 

Requirements  3 and  4 lead  some  people,  in- 
cluding the  author,  to  view  the  logarithmic  profile 
as  applying  basically  to  area-average  profiles, 
rather  than  to  point  profiles.  Requirements  3 
and  4 are  then  looked  upon  as  requirements 
for  obtaining  area-average  profiles  with  a minimum 
of  instrumentation. 

The  interpretation  put  upon  the  three  param- 
eters— F* , Zo  (sometimes  called  roughness  length) , 
and  D (sometimes  called  displacement  height) — 
has  been  extremely  diverse  during  the  past  24 
years.  Extreme  views  have  been  that  roughness 
length  is  “a  constant  of  integration’’  on  the  one 
hand,  and  a physical  parameter  of  the  ground 
surface  on  the  other.  The  parameter  F*  has 
generally  been  identified  as: 


Recent  work  shows  that  for  flexible  elastic  plant 
covers  such  as  corn  or  rice,  both  Zq  and  D vary 
with  windspeed,  as  indicated  by  F,^,  geostrophic 
wind,  or  speed  at  a reference  height.  The  study 
of  the  manner  of  this  dependence  would  be  facili- 
tated by  a convenient  means  of  determining  F,^, 
Zq,  D,  and  their  standard  errors  from  observations. 

On  receipt  of  some  notes  from  Stephen  M. 
Robinson  and  C.  B.  Tanner  ^ at  the  University 
of  Wisconsin,  giving  substantially  the  material 
in  pages  1 to  5 of  Technical  Note  No.  1 {26), 
a program  for  machine  computation  of  adiabatic 
wind  profile  parameters  was  written  here.  This 
program  accomplishes  the  same  data  reduction 
task  as  Mr.  Robinson’s,  and  is  based  on  equations 
written  by  him.  It  has  two  major  differences 
in  technique  that  may  prove  useful.  The 
subprogram  for  standard  errors  and  the  equa- 
tions on  which  they  are  based  are  original  con- 
tributions. The  program  is  written  in  the  Bur- 
roughs version  of  Algol,  for  the  Burroughs  220 
computer.  This  report  discusses  mathematical 
methods  and  not  the  details  of  the  computer 
program. 

REVIEW  OF  THE  PROBLEM 

The  equation  to  be  fitted  to  the  adiabatic 
wind  profile  observations  is: 


V.=iVr/p 


(86) 


where  k is  Von  Karman’s  constant,  ~0.4;  r 
is  shear  stress  in  the  surface  layer;  and  p is  air 
density.  Let  Zs  be  the  value  of  2:  at  which  equa- 
tion 85  gives  F=0.  Then: 


or. 


0=F*  In 


z,=  -D^Zo.  (87) 


The  parameter,  D,  is  therefore  related  to  the 
height  at  which  the  profile  extrapolates  to  zero, 
differing  by  Zq,  and  having  a change  in  sign. 

For  rough,  rigid  surfaces,  Zq  appears  to  be  a 
physical  parameter  of  the  surface  and  2^  = 0. 
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Robinson’s  development  is  followed  in  general, 
with  summation  and  mean  value  notation  rather 
than  vector  notation.  Data  are  speeds  Ff  for 
the  n measuring  heights,  Zi. 

The  sum  of  the  squares  of  errors  to  be  minimized 
in  evaluating  F,^,  D,  and  Zq  is: 

i?=±(F<-nin(£l±^)J.  (88) 

Let: 

W=ln  Zq, 

Xi=\n  {Zi+B), 


^ Personal  communication. 
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yi=Xi-x, 

u,=V,-V, 

ri=ll{Zi-\-D). 

Put  the  three  partial  derivatives  simultaneously 
equal  to  zero : 

^^=-22  {[V,-V,(x,-w)]{x,-w)]=0-,  (89) 

OV  ^ 1=1 

A 77"  ^ 

^=-2V,T.{[Vi-V^{,x,-w)]r,]=Q-,  (90) 

OU  i = i 

^=‘^±,{V,-V^(x,-w)]=0.  (91) 

^ OZq  Zq  2=1 

In  equations  89,  90,  and  91,  w is  a function  of  2o 
and  the  Xi  and  depend  upon  D.  The  set  of 
values  (F*,  D,  Zo)  is  required  that  satisfies 
equations  89,  90,  and  91. 

In  equation  91,  averages  are  taken,  and  the 
equation  is  rearranged  to  give: 


A form  of  the  latter  function  is  chosen  by 
Robinson.  The  choice, 


_ F 

W=X  — ^' 
V * 


a(Z))=constant  • Xi— 

\y2  uy-yrf 


(97) 


(92) 


has  the  merit  that  g{D)  so  defined  has,  for  reason- 
able data,  a nearly  linear  form  for  a wide  range  of 
trial  values  of  D.  This  near-linearity  is  convenient 
when  new  estimates  of  the  root,  D,  are  found  by 
linear  interpolations.  Hence,  equation  97  is  used. 

The  second  difference  in  technique  is  in  finding 
the  first  two  and  subsequent  estimates  of  D.  In- 
stead of  a pair  of  reasonable  guesses  being  supplied, 
the  program  automatically  makes  two  extreme 
guesses,  to  bracket  the  root.  New  estimates  are 
found  alternately  by  linear  interpolation  and  by 
bisection,  until  one  of  three  criteria  calls  for  making 
a final  linear  interpolation.  Each  new  estimate  is 
used  to  narrow  the  bracket. 

The  data  format  consists  of  integers  (groups  of 
digits)  separated  by  spaces: 

aaa  bb  ccc  ccc  ccc  ccc. 

bb  integers 


Equations  90  and  89  are  rewritten  and  averaged, 
and  the  value  of  w from  equation  92  is  substituted 
in  them.  This  gives: 

{)=W—V^  Ijr  (93) 

and  

0=wy—V^  y\  (94) 

From  equation  93  comes  the  equation  for  F,,; : 

F,==-  (95) 

yr 

Combine  equations  94  and  95  to  obtain: 

0=WJ-^-  (96) 

yr 

Equation  96  is  the  implicit  equation  in  the  single 
unknown,  D,  which  is  to  be  solved.  This  is  done  by 
defining  a function,  g{D),  which  has  as  a single 
real  root  the  desired  value  of  D,  and  by  finding 
that  root.  Any  of  several  functions  g{D)  will  do: 

(7'(/7)=constant  • {uy—y^  wlyr] 
or 

^"(75)  = constant  • [y^  ur—uy  Ipr] , 
for  example. 


If  the  first  integer  (aaa)  then  (aaa)  is  the 
profile  identification  number,  and  the  (ccc)  are 
velocities,  starting  at  the  lowest  anemometer 
height.  The  second  integer  (bb)  is  the  number  of 
measuring  heights.  If  (aaa)=0,  then  the  (ccc) 
are  anemometer  heights,  beginning  with  the  lowest. 
The  computer  subsequently  converts  input  data  to 
floating  point  form,  and  the  output  has  explicit 
decimal  points.  Units  of  F,^  in  the  output  are 
those  of  Vi  in  input;  units  of  D and  Zq  in  the  output 
are  those  of  Zi  in  input.  Units  of  centimeters  per 
second  and  centimeters  were  visualized  when 
writing  the  program,  but  are  not  essential. 

The  first  and  second  trial  values  of  D are 
Z>i  = — 2i  + l and  P2=  + 2i  — 1.  If  these  do  not 
establish  a bracket  for  D,  the  fact  is  noted  and  the 
next  profile  is  read  in. 

The  three  criteria  for  making  a final  estimate  of 
D by  linear  interpolation  are: 

(1)  The  absolute  value  of  g{D)  less  than  a 

critical  amount,  or 

(2)  The  size  of  the  bracket  less  than  a critical 

amount,  or 

(3)  The  number  of  trial  values  of  D reaches  a 

maximum  allowed  number. 

With  several  sets  of  real  data,  all  profiles  having 
four  or  more  measuring  heights,  and  heights  being 
measured  from  approximately  the  mean  soil  sur- 
face, seven  trial  values  of  D (including  the  first 
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two)  were  always  sufficient  to  evaluate  the  dis- 
placement height,  D,  to  the  nearest  hundredth  of 
a centimeter,  if  the  least-squares  value  was  within 
the  original  bracket  (fig.  18). 

A SAMPLE  COMPUTATION 

The  same  data  are  used  as  in  Robinson’s  sample 
calculations.  These  are  Lettau’s  mean  adiabatic 
profile  from  the  Johns  Hopkins  observations  at 
O’Neill,  Nebr.,  in  1953  {15). 


2 ( cm.)  V ( mm. ! sec.) 

40 4,902 

80 5,853 

160 6,737 

320 7,592 

640 8,375 


Successive  estimates  of  D and  corresponding 
v^alues  of — 


!7(0)=i 

y" 


ur 

^yr 


are: 


Method:  D g{D)  Resulting  bracket 

-39.0  +0.029518  

+ 39.0  -0.060772  -39.0  +39.0 

By  linear  interpolation —13.499818  +0.006182  —13.5  +39.0 

By  bisection +12.750091  -0.030611  -13.5  +12.75 

By  linear  interpolation —9.088715  —0.000678  —13.5  —9.09 

By  bisection -11.294266  +0.002713  -11.29  -9.09 

By  linear  interpolation —9.  529773  —0.  000006  —11.  29  —9.  529773 


Then  a final  linear  interpolation  gave  —9.53  cm.  approximately  linear  dependence  of  g{D)  on  D for 

for  D.  The  corresponding  Zq  was  0.432  cm.,  and  the  example  above. 

V:,,  was  1,151  mm. /sec.  Figure  18  shows  the 


ESTIMATING  THE  ACCURACY  OF  V^, 
D,  AND  Zo  FROM  ADIABATIC  WIND 
PROFILE  DATA 


^Vi  I v^,  Vs ^»-  H 


dUi 


Vo.  Vo V„.  B 


m 

dD 

Vi,  V2.  . 

• • 

(100) 


Using  a computer  with  eight  significant-figure 
arithmetic  does  not  lead  to  computational  error 
in  the  profile  parameters.  Errors  depend  upon 
the  data  for  the  heights,  Zi,  and  the  speeds,  V^.  The 
heights  can  be  assumed  to  be  accurate,  and  all  the 
error  of  measurement  (sampling  and  instrumental) 
can  be  assigned  to  the  speeds.  Then  statistical 
analysis  leads  to  estimates  of  the  errors  in  com- 
puted parameters  U,,,  D,  and  Zq  as  functions  of 
the  data.  This  statistical  analysis  involves  only 
the  assumptions  already  implied  in  the  fitting  of  a 
least-squares  profile. 

The  Error  in  D 

Define 


and,  of  course,  similar  equations  apply  with  respect 
to  V2,  Vs,  . . . Vn.  The  error  in  displacement 
height,  D,  is  given  by: 

E(y.)+ . . EiV„) 

+ (higher  order  terms) 

where  the  derivatives  are  evaluated  at  the  true 
values  of  the  velocities. 

The  standard  error  of  D is  given  (approximately) 

= (101) 


H{D,  Vi,  V2,  . . . Vn)=u7/  yr  — ip  ur  (98) 

where  D is  an  arbitrary  variable,  and  the  Vi  are 
measured  velocities. 

Then, 


If  the  standard  error  of  U<  is  assumed  constant 
with  height,  then: 

dD 


dH — dD-\-^2i  vfT' 

dD  y Vo V,  i^idVi 


dVi.  (99) 


Equations  are  required  for 


Divide  equation  99  by  dVi,  and  put  0=dV2= 
dVs=  . . . dVn=dB.  Then: 


Vj,  V2, 


and 

v„ 
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Figure  18. — The  functions  g{D)  for  the  sample  computation. 

SO  that 


H=uy  yr—y^  ur, 


therefore: 

dH  duy  — , dyr  dy^ dur 


or 


yr-\-uy  sr—uy  yr^—2  ys  ur-\-y^  ur^, 

(103) 

where  Si=r,—r.  Evaluating  at  the  true  velocities 



SO  that Ui=V^yi,  gives  ■^=V^{y^  sr—ys  yr). 


Similarly, 


dH  yt  -5  Si 

-^  = yr  ^ — -1, 

dFv  n ^ n 


dP 

dVi 


r Si—yr  yi 


Vi9^v  H 


dH 

dP 

v„  v„  . 

■ • 

The  Error  in  F* 


ur 


is  computed  from 
The  standard  error  of  is: 


(lOo) 


But, 


Vi9^i  ^P 


dP  I dF, 

Vi  dVj  lvi^i  dVj 


(106) 
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—ur^  yr — ur  rs + ur  bD 

(P)2  dVj 


+4l. 

,•  nyr 
(107) 


Evaluating  at  the  true  value  of  velocities. 


dF, 


rs  bD 
yr  bV j 


nyr 


The  Error  in  Zq 

The  roughness  length,  Zq,  is  given  by: 
2„=exp  (z-^y 

he  standard  error  of  Zq  is  approximated 


But, 


^^0 

dF, 


This  leads  to 


^^0 

bD 


bP 

Vi  bV j 


I ^^0 

aF, 


(108) 


(109) 


(110) 


Wi 


■ur  rs 


VbF,k^,.^F*L  njj 

If  we  evaluate  at  true  values. 


-\-ur  yr^ 


(111) 


bzQ 

bF, 


/ /_  V—s\  ^1 


nV^yr  n)} 


To  compute  these  standard  errors  of  the  profile 
parameters,  several  intermediate  variables  are 


required  in  addition  to  those  found  in  computing 
the  profile  parameters.  These  are^for  tlm  least 
squares  values  pf  D,  F,^:,  and  Zq\  r,  'sr,ys,T2;  and, 
at  each  measuring  height,  i:  Si  =ri  — r; 


bP 

bVi 


H,  Vj^i 


from  equation  100; 


from  equation  107; 


from  equation  111. 


bF, 


bgp 

bF, 


Viy^i 


Date 

10  July 

11  July 

23  July 

24  July 

24  July 

Time 

1905 

0605 

1905 

0605 

0705 

Tl6m  — T.25in 

+ 0.24°C 

+ 0.  23 

-0.  07 

+ 0.  48 

-0.  92 

Ul6m 

729  cm/sec 

780  cm/sec 

731  cm/sec 

584  cm/sec 

843  cm/sec 

Ugm 

639 

689 

635 

532 

763 

U4m 

569 

607 

568 

505 

708 

U2m 

497 

543 

420 

609 

Ulm 

432 

474 

437 

380 

552 

U.5m 

373 

411 

377 

336 

501 

U.25m 

299 

340 

310 

289 

426 

V*,  cm/sec 

108  ±4 

113±5 

108  ±6 

75±6 

107  ±6 

Zo,  cm. 

2.  0±0.  5 

1.  7±0.  4 

2.  0±0.  7 

. 65±.  38 

. 61±.  26 

D,  cm 

9.  0±5.  3 

12.  0±6.  0 

12.  4±8.  5 

5.  3 ±10.  3 

9.  6±8.  0 

The  standard  error  of  the  windspeeds  is  also 
required,  SE{Vy).  This  may  be  estimated,  for 
n>4,  by 

■ {u^-V^wyY'K  (112) 

In  any  case,  SE{Vi)  is  at  least  as  large  as 
known  anemometer  errors.  The  parameters  for 
the  sample  problem  are  now  repeated,  with  their 
standard  errors  indicated: 

F^  = 115.1±1.4  cm./sec. 

2Q=0.43dr0.04  cm. 

P—  —9.5±1.2  cm. 

ANALYSIS  OF  A SET  OF  NEAR- 
ADIABATIC  PROFILES 

For  a set  of  good  adiabatic  wind  profiles,  six 
measured  at  O’Neill,  Nebr.,  in  the  summer  of 
1956  by  Halstead  and  associates  (6)  were  chosen. 
Data  and  results  of  analysis  are  tabulated  below. 


24  July 
1905 
-0.  18 
503  cm/sec 
427 
398 
339 
291 
259 
224 


>24 
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Compared  to  the  mean  profile  first  discussed, 
these  six  have  appreciably  higher  standard  errors 
for  the  parameters.  That  D-Zq  is  positive  is 
a minor  puzzle.  It  indicates  that  the  windspeeds 
at  the  highest  levels  are  less  than  expected,  while 
those  at  the  lowest  levels  are  greater.  At  the 
highest  level,  the  log-law  for  area  averages  does 
not  hold  so  well;  at  the  lowest  levels,  the  measure- 
ments are  not  sufficiently  representative  of  area 
averages.  In  this  instance,  the  tendency  to  avoid 
local  obstacles  and  depressions  in  locating  the 
anemometer  mast  appears  to  have  led  to  selecting 
a spot  with  slightly  higher  winds  than  average 
near  the  surface.  (These  data  were  not  taken 
with  this  manner  of  statistical  profile  analysis  in 
mind.  Subjective  analysis  of  profile  data  gen- 
erally weights  the  anemometer  readings  according 
to  expected  reliability.)  These  new  computations 
show  again  that: 

(a)  Site  selection  and  anemometry  were  good ; 


(b)  The  magnitude  of  displacement  height, 
D,  was  small  compared  to  median  height 
of  measurement;  and 

(c)  Zo  was  approximately  1 cm. 

SUMMARY 

Another  program,  based  on  Robinson’s  equa- 
tions (^6),  has  been  written  for  evaluating 
logarithmic  wind  profile  parameters.  Two  fea- 
tures in  the  technique  differ  from  Robinson’s:  (1) 
the  method  of  choosing  the  first  estimates  of  D, 
and  (2)  the  method  of  finding  subsequent  estimates 
of  D.  These  perhaps  are  improvements.  In 
addition,  this  program  computes  good  estimates 
of  the  standard  errors  of  the  three  parameters, 
provided — 

(a)  The  accuracy  of  the  anemometers  is 
Imown,  or 

(R)  A sufficient  number  of  measuring  heights 
are  used. 


TURBULENT  TRANSFER  CHARACTERISTICS  OF  THE  AIRSTREAM 
IN  AND  ABOVE  THE  VEGETATIVE  CANOPIES  AT  THE  EARTH’S 

Surface 


By  Jerry  Stoller  and  Edgar  R.  Lemon 


Turbulent  exchange  calculations  for  sensible 
heat,  latent  heat  (evaporation  and  condensation), 
and  carbon  dioxide  (photosynthesis  and  respira- 
tion) are  important  in  estimating  some  of  the 
components  of  the  energy  balance  near  the  earth’s 
surface.  Basic  to  these  calculations  is  the  deter- 
mination of  the  transfer  coefficient  for  momentum 
{Km)  and  applying  it  to  the  transfer  of  other 
molecular  properties.  In  essence,  one  implicitly 
determines  the  value  of  Km  in  the  turbulent 
boundary  layer  above  the  surface  and  uses  it 
to  calculate  the  exchange  of  other  molecular 
properties,  knowing  their  potential  drop  (gradient) 
across  some  known  distance  in  the  same  boundary 
layer.  The  assumption  is  made  that  the  momen- 
tum transfer  coefficient  is  quantitatively  equal 
to  the  other  transfer  coefficients,  i.e.,  for  heat, 
water  vapor,  and  carbon  dioxide. 

Upon  the  accurate  measurement  of  the  momen- 
tum transfer  characteristics  of  the  turbulent 
boundary  layer,  then,  rests  the  whole  aerodynamic 
approach  to  the  energy  budget. 

TURBULENT  EXCHANGE  ABOVE  THE 
CROP  CANOPY 

Where  vegetation  exists  at  the  earth’s  surface, 
the  turbulent  boundary  layer  has  to  be  sub- 
divided into  two  regions:  One  above  the  plant 
surfaces  where  no  sources  or  sinks  occur  and  the 
other  within  the  vegetative  canopy  where  sinks 
and  sources  do  occur.  Exchange  calculations 
above  the  canopy  rely  upon  the  log  wind  profile 
law  and  the  logarithmic  law  of  vertical  distribu- 
tion of  the  other  physical  quantities.  Three 
specific  problems  arise,  however,  in  the  use  of 
this  aerodynamic  approach  to  exchange  calcula- 
tions above  the  vegetative  canopy: 

(1)  The  errors  created  by  appreciable  thermal 
gradients  in  the  turbulent  boundary  layer; 

(2)  Instrument  and  interpretative  limitations 
at  low  mean  wind  velocities;  and 

(3)  The  functional  relationships  of  the  log 
wind  profile  characteristics  to  both  mean  wind- 
speed  and  surface  properties  (geometric  and  elas- 
tic properties  of  the  roughness  element) . 

When  the  temperature  of  the  air  near  the 
vegetation  rises  significantly  during  the  daytime, 
then  thermal  convective  transfer  processes  (“free 
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convection”)  become  important  in  addition  to  the 
mechanical  convective  transfer  processes  (“forced 
convection”).  [The  forced  convection  processes 
are  responsible  for  the  log  profile  relationship.] 
Free  convection  superimposed  on  forced  convec- 
tion causes  an  underestimation  of  the  transfer 
processes  when  the  log  profile  law  is  used,  and  is 
especially  pronounced  at  low  wind  velocities  when 
the  forced  convection  component  is  small.  Some- 
times the  mean  wind  velocity  is  reduced  to  such 
a low  level  that  it  is  not  measurable  by  standard 
means  (rotating  cup  anemometers).  Under  these 
conditions,  the  logarithmic  law  of  vertical  distri- 
bution of  physical  quantities  no  longer  holds 
true  and  the  calculated  fluxes  erroneously  go  to 
zero.  The  subject  of  “free  convection”  error 
will  be  dealt  with  in  detail  in  a later  report. 

In  addition  to  free  convection  and  low  wind- 
speed  problems  associated  with  the  log  law  ap- 
plication to  transfer  processes  in  the  turbulent 
airstream,  another  problem  arises  out  of  the  func- 
tional relationship  between  ^vindspeed  and  the 
wind  profile  characteristics  over  vegetative  sur- 
faces. This  problem  is  created  by  the  changing 
geometric  properties  of  the  surface  (plant  growth) 
and  the  changing  elastic  properties  of  the  surface 
(waving  of  stalks  and  leaf  flutter).  This  func- 
tional relationship  will  be  taken  up  in  detail 
first. 

The  eddy  flux  across  a plane  of  some  entity 
(i.e.,  momentum)  which  is  neither  destroyed  nor 
created  in  its  transfer  (not  strictly  true  for  momen- 
tum) can  be  written: 

(113) 

where  Q is  the  total  amount  of  the  entity  passing 
downward  through  a unit  area  per  unit  time ; K is 
the  transfer  coefficient,  and  bcj^z  is  the  gradient 
across  the  plane.  The  amount  of  the  entity,  c, 
is  expressed  on  a per -unit-volume  basis,  and  2 is 
the  axial  distance  normal  to  the  plane. 

For  the  flux  of  momentum  (shearing  stress), 
T,  the  basic  equation  becomes: 

T=Km^p  (114) 

where  u is  the  windspeed  and  p is  the  density  of 
the  air. 
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In  the  turbulent  boundary  layer  immediately 
above  the  vegetation  surface,  the  vertical  gradient 
of  wind  is  expressed  by : 


bz  k\p  ) z 


(115) 


and  upon  integration  gives  the  “logarithmic  wind 
profile  law”  where  the  windspeed,  at  height 
0 above  the  ground  is  represented  by: 

^ In  where  z'^Zq,  (116) 

To  is  the  shearing  stress  at  the  surface,  Zq  is  the 
roughness  length,  and  k is  the  Von  Karman 
constant  (0.4). 

Equations  115  and  116  imply  that  the  windspeed 
profiles  are  functions  of  log  height,  s,  above  a 
given  reference  plane.  As  the  ground  vegetation 
grows,  however,  the  turbulent  boundary  layer  is 
displaced  upward,  causing  the  reference  plane 
from  which  s is  measured  to  be  displaced  upward 
also  from  the  ground  surface.  It  then  becomes 
necessary  to  introduce  a zero  point  displacement, 
d,  and  rewrite  equation  116,  thus: 


gradients  were  small).  Three  conditions  mini- 
mizing this  error  were  observed:  (1)  The  measure- 
ments were  made  near  the  surface  (within  160  cm. 
of  the  ground,  or  85  cm.  of  the  crop  top) ; (2)  the 
measurements  were  made  during  cloudy  periods 
when  the  radiation  load  was  small;  and  (3)  the 
alfalfa  was  lush  with  plenty  of  soil  moisture 
(soil  moisture  tension  was  0.2  atm.  at  6-inch 
depth) . 

Table  1 and  figure  19  present  the  functional 
relationships  between  the  mean  windspeed  at 
160  cm.  and  the  profile  parameters  Zo,  D,  and  V*. 
It  can  be  seen  that  the  roughness  length,  Zq,  ranged 
from  28.5  to  1.3  cm.,  and  the  effective  displace- 
ment, D,  from  —18.8  to  —65.0  cm.  The  friction 
velocity,  V*,  appeared  fairly  independent  of  wind- 
speed.  (Note  that  table  1 and  figure  19  give 
values  of  V*/k  instead  of  V*.)  This  would  indi- 
cate that  the  shearing  stress  was  fairly  constant 
and  independent  of  the  windspeed  at  a given 
height  above  the  alfalfa.  One  can  conclude  that 
within  the  range  of  windspeeds  encountered, 
shearing  stress  is  almost  completely  dependent 
upon  the  physical  characteristics  of  the  surface. 


+ 

2o 


(117) 


where  D=d-]-ZQ,  or  the  “effective  displacement 
parameter,”  and  Zi  is  the  nominal  height  above 
the  ground  surface. 

The  friction  velocity,  V*,  defined  by: 


then: 


and, 


ku 


17*=- 


Zi-\-D 


In 


^0 


Km  = {zi-^D)kV^. 


(118) 

(119) 

(120) 


The  preceding  section  (p.  28)  takes  up  the 
method  used  here  for  the  determination  of  the 
logarithmic  wind  profile  parameters,  V*,  Zq,  and 
D from  sufidcient  experimental  data  u on  zt. 

In  that  section,  V*. 

k 

Alfalfa  Investigations 

Numerous  wind  profile  measurements  over 
alfalfa  were  taken  with  Sonoya  Rotating  Cup 
Anemometers  on  July  29,  30,  and  31,  1960,  when 
the  alfalfa  was  75  cm.  high.  The  equipment  and 
site  have  been  fuUy  described  elsewhere  (fj). 
Mean  wind  velocities  were  measured  over  5- 
minute  periods  at  90,  100,  120,  and  160  cm.  from 
the  ground.  Sampling  times  were  selected  when 
it  was  felt  that  thermal  effects  upon  the  profiles 
would  be  negligible  (i.e.,  when  temperature 
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Figure  19. 


-Wind  profile  parameters  for  alfalfa,  July 
29-31,  1960. 
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Table  1. — Functional  relationships  between  windspeed  and  roughness  length  (so),  effective  displacement 
(Z>),  and  friction  velocity  (V*)  for  alfalfa  at  Ithaca,  N.Y.,  July  29-31,  1960 


1 

D 

V*lk 

n. 

15.  6 

Cm. 

-48.  4 

Cm. /sec. 

88.  8 

25.  2 

-23.  3 

121.  2 

28.  5 

-18.  8 

140.  9 

17.  8 

-32.  8 

118.  2 

16.  4 

-38.  8 

122.  8 

11.  0 

-46.  4 

114.  7 

11.  3 

-39.  8 

117.  6 

3.  5 

-57.  8 

94.  6 

10.  7 

-33.  5 

141.  7 

3.  2 

-59.  0 

no.  4 

3.  2 

-57.  9 

118.  1 

5.  3 

-49.  4 

137.  6 

2.  8 

-57.  8 

125.  0 

1.  3 

-65.  0 

106.  2 

2.  7 

-55.  8 

137.  1 

Profile  No. 


Date  and 
time 


Height 


Velocity 


1_ 

2_ 

3_ 

4_ 

5_ 

6_ 

7_ 

8_ 

9. 

10 

11 

12 

13 

14 

15 


E.s.t. 

7/30/60 

1027 


7/30/60 

1002 


7/30/60 

0956 


7/30/60 

1206 

7/29/60 

1840 

1852 

1904 

7/29/60 

1916 

7/30/60 

0950 

7/30/60 

1008 

1014 


7/29/60 

1846 


7/29/60 

1605 


7/31/60 

1259 


7/31/60 

1311 

7/31/60 

1231' 

1353 


7/31/60 

1317 


7/31/60 

1247 

7/31/60 

1219) 

12531 


Cm. 


90 

100 

120 

160 

90 

100 

120 

160 

90 

100 

120 

160 

90 

100 

120 

160 

90 

100 

120 

160 

90 

100 

120 

160 

90 

100 

120 

160 

90 

100 

120 

160 

90 

100 

120 

160 

90 

100 

120 

160 

90 

100 

120 

160 

90 

100 

120 

160 

90 

100 

120 

160 

90 

100 

120 

160 

90 

100 

120 

160 


Cm.  I sec. 

87 

107 

135 

175 

118 

135 

163 

205 

130 

145 

180 

225 

138 

156 

188 

232 

140 

160 

197 

245 

158 

182 

218 

268 

175 

197 

230 

278 

212 

235 

275 

320 

235 

260 

295 

350 

250 

280 

325 

380 

270 

305 

347 

408 

280 

315 

355 

420 

305 

340 

387 

450 

312 

350 

395 

455 

345 

385 

430 

500 
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When  one  observes  the  motion  of  the  alfalfa 
plants  in  the  field,  two  factors  readily  become 
apparent.  As  the  wind  velocity  increases,  the 
leaves  orient  themselves  in  the  direction  of  the 
wind,  and  the  field  appears  to  consist  of  waves; 
i.e.,  “smoothing  out  effect.” 

Under  low  velocities,  the  plants  stand  nearly 
erect  while  each  leaf  acts  as  a momentum  sink. 
The  momentum  loss  is  probably  greatest  in  the 
area  of  the  upper  one-third  of  the  plants;  however, 
some  wind  movement  is  apparent  down  to  the 
ground  level.  As  wind  velocities  increase,  how- 
ever, the  upper  leaves  orient  themselves  in  the 
direction  of  the  wind  flow  and  the  plants  overlap 
each  other.  Thus,  a complete  surface  of  smooth 
leaves  is  approached  and  is  accompanied  by  a drop 
in  the  roughness  length,  Zq-  For  the  effect  to 
occur,  the  plant  must  be  high  enough  to  provide 
a complete  ground  cover  and  consist  of  enough 
foliage  to  create  the  “wave”  appearance  under 
higher  velocities. 

Wheat  Investigations 

The  air  turbulence  studies  in  a wheatfield  were 
made  on  an  open,  nearly  level  hilltop  (Mount 
Pleasant)  located  adjacent  and  to  the  north  of 
Ellis  Hollow  Experimental  site  described  else- 
where (14)-  The  field  has  the  approximate 
dimensions  of  200  by  200  meters.  As  in  the 
alfalfa  investigations,  measurements  were  made 
under  conditions  favoring  small  temperature 
gradients,  and  with  the  same  equipment. 

It  was  observed  earlier  that  within  the  custom- 
ary 5-minute  averaging  period  there  were  periods 
of  relatively  steady  wind.  In  an  effort  to  separate 
the  steady  flow  from  the  gusts,  a series  of  30- 
second  sampling  periods  were  used.  When  a 
series  of  two  or  more  consecutive  30-second  periods 
rendered  similar  readings  from  the  cup  anemom- 
eters, the  profile  was  considered  to  be  semi-steady. 
By  selecting  semi-steady  periods  the  semi-log 
plotting  of  the  wind  profiles  proved  to  be  much 
more  precise  than  when  the  5-minute  averaging- 
period  was  employed. 

On  the  first  day  of  study  (July  9,  1960),  the 
wheat  averaged  almost  100  cm.  in  height  at  a 
density  of  60  stalks  per  square  foot.  The  heads 
were  10  cm.  long,  and  located  right  above  the 
upper  leaves.  Wind  profiles  were  taken  when  the 
wind  blew  from  a northeasterly  direction.  As  can 
be  seen  in  figure  20  and  table  2,  V*  and  Zq  in- 
creased with  an  increase  in  wind  velocity.  Also, 
the  effective  displacement,  D,  decreased  when  the 
wind  velocity  above  the  crop  increased. 

In  regard  to  these  measurements,  it  must  be 
rnentioned  that  the  wind  velocities  did  not  attain 
high  enough  values  to  cause  appreciable  bending 
of  the  wheat  stalks.  A different  picture  was 
obtained  July  22  and  23,  when  measurements 
were  taken  at  the  same  site.  At  this  time,  the 


average  height  of  the  wheat  was  approximately 
130  cm.  Most  of  the  growth  had  taken  place 
between  the  upper  leaf  and  the  head.  Thus,  the 
center  of  mass  was  displaced  upward. 

Figure  21  and  tables  3 and  4 present  the  data 
showing  that  the  roughness  length,  Zq,  and  effec- 
tive displacement,  D,  decreased  while  the  friction 
velocity,  U*,  remained  fairly  constant  with  in- 
creasing windspeed  above  a “critical”  level.  Evi- 
dently, the  wheat  crop  on  July  22,  1960,  acted 
similarly  to  alfalfa.  The  reversal  in  trends  for  the 
wheat  crop  during  the  two  different  stages  of  de- 
velopment has  to  be  attributed  both  to  a change- 
in  the  physical  characters  of  the  crop,  and  to  the 
diversity  of  range  of  windspeeds  encountered.  In 
the  earlier  stage  (table  2),  the  winds  were  slower 
and  the  crop  shorter,  thus  more  rigid.  This  would 
permit  some  waving  and  leaf  flutter,  but  little 
streamlining.  At  the  later  dates  under  the  higher 
wind  velocities  when  the  crop  was  taller,  stream- 
lining probably  occurred  much  as  it  did  in  alfalfa. 
It  is  of  interest  to  point  out  that  the  trends  ob- 
served in  the  wheat  at  the  lower  windspeeds  on 


c m./sec. 


«170,„,  (cm./sec.) 

Figure  20. — Wind  profile  parameters  for  wheat,  July  9, 
1960. 


38 


PRODUCTION  RESEARCH  RPT.  72,  U.S.  DEPT.  OF  AGRICULTURE 


Table  2. — Wind  'profile  data  jor  'wheat  at  Ithaca, 
N.Y.,  July  9,  1960 


Profile 

No. 

Height 

Velocity 

V*lk 

D 

2c 

Cm. 

Cm.jsec. 

Cm.jsec. 

Cm. 

Cm. 

( 100 

180 

\ 

110 

225 

1 

1 130 

310 

> 125.  3 

-82.  8 

4.  20 

170 

380 

1 250 

460 

J 

( 100 

250 

110 

345 

2 

1 130 

420 

> 108.  1 

-91.  9 

0.  788 

170 

485 

[ 250 

580 

) 

f 100 

320 

\ 

no 

390 

3 

{ 130 

455 

> 195.  6 

-68.  0 

6.  03 

170 

545 

1 250 

670  , 

) 

( 100 

400 

\ 

no 

470 

4 

< 130 

535 

) 266.  5 

-49.  4 

10.  9 

170 

630 

1 250 

780  , 

) 

(cm. /sec.) 


300  400  500  600  700  800  900 

~170cm.(cm./sec.) 


the  earlier  date  parallel  those  found  in  mature 
corn  at  higher  windspeed  in  our  1959  investiga- 
tions {14)-  A few  values  for  corn  of  F*,  Zq,  and 
as  a function  of  windspeed  are  given  later  in 
table  5.  They  are  in  agreement  with  the  1959 
trends. 

Table  3. — Wind  profile  data  Jor  'wheat  at  Ithaca, 
N.Y.,  July  22,  1960 


1 


2 


3 


4 


5 


Profile 

No. 


Height 


Velocity 


V*/k 


D 


20 


Cm. 


Cm.jsec. 


Cm.jsec. 


Cm. 


Cm. 


[ 140 

150 
{ 170 

210 
I 290 
( 140 

150 

< 170 
210 

I 290 
( 140 

150 

< 170 
210 

I 290 
( 140 

150 

< 170 
210 

I 290 
( 140 

150 
I 170 
210 
I 290 


380 

410 

440 

500 

585 

440 

470 

500 

585 

670 

550 

585 

640 

735 

840 

585 

640 

725 

840 

990 

640 

725 

820 

950 

1080 


I 210.  2 
I 221.  6 

I 228.  7 
I 285.  8 
I 213.  0 


-46.  6 


-58.  8 


-82.  2 


-91.  6 


-118.  0 


15.  0 


11.  1 


5.  24 


6.  22 


1.  09 


Figure  22  gives  a qualitative  diagram  of  the 
functional  relationships  for  the  three  crops  for 
comparative  purposes.  The  trends  on  wind- 
speed  for  the  different  crops  are  chiefly  a matter  of 
rigidity  of  momentum  sink  elements,  whereas 
the  order  of  magnitude  of  the  values  is  attributed 
to  the  momentum  sink  geometry.  The  alfalfa 
is  the  least  rigid,  the  taller  wheat  next,  followed 
by  the  shorter  wheat,  with  corn  being  the  most 
rigid.  A maximum  Zq  probably  takes  place 
in  all  three  crops.  Low  enough  windspeeds  in 
alfalfa  were  not  encountered  to  clearly  define 
this  point,  whereas  it  might  take  very  high  winds 
to  observe  a reversal  in  corn.  The  wheat  would 
fall  between  the  alfalfa  and  corn  in  characteristics. 
Of  course,  all  these  relationships  would  depend 
upon  stage  of  crop  maturity,  plant  density,  and  a 
host  of  other  factors  too  numerous  to  mention. 

TURBULENT  EXCHANGE  WITHIN  THE 
PLANT  CANOPY 


Figure  21. — Wind  profile  parameters  for  wheat,  July 
22-23,  1960. 


The  exchange  processes  within  the  vegetative 
canopy  are  theoretically  much  more  complex  than 
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(cm.) 


FiGrRE  22. — Idealized  roughness  parameters  for  alfalfa,  AA^heat,  and  corn  as  a function  of  windspeed. 


those  above,  because  of  the  complex  distribution  of 
sources  and  sinks  within  the  “crop  volume.” 
Instrument  and  interpretative  limitations  on  the 
experimental  level  are  restrictive,  too.  None- 
theless, we  have  attempted  to  arrive  at  a semi- 
quantitative  picture,  at  least,  of  some  transfer 
characteristics  within  the  wheat  and  corn  canopies. 

Referring  again  to  equation  113,  we  have 


which  states  that  the  downward  flux,  Q,  of  some 
entity  is  equal  to  the  product  of  the  transfer  co- 
efficient, K,  and  the  gradient  of  the  entity  across 
some  plane,  dcidz.  If  we  assume  that  the  transfer 
coefficient  for  momentum  is  quantitatively  equal 
to  the  transfer  coefficient  for,  say  CO2,  water 
vapor,  and  heat,  then  experimentally  we  need  to 
determine  the  transfer  coefficient  for  momentum 
from  windspeed  characteristics  and  the  gradient 
for  the  entity  in  question  (i.e.,  CO2,  water  vapor 
or  heat) , in  order  to  calculate  the  exchange  rates 
for  the  latter. 

Using  Prandtl’s  model  for  rough  flow  in  a pipe 
(33)  and  making  the  assumptions  that  turbulence 
is  fully  developed  and  isotropic  within  the  vege- 
tative canopy,  we  have: 


(121) 


where  K(Zi)  is  the  momentum  transfer  coefficient 
at  nominal  height,  2^;  I is  the  “mixing  length”; 
and  w'  is  the  vertical  windspeed  fluctuation,  or 
“vertical  eddy  velocity.”  The  eddy  velocity,  w', 
is  defined  as  the  difference  between  the  instan- 
taneous velocity,  w,  and  a mean  velocity,  w. 

Experimentally,  w'  and  I were  evaluated  from 
hot  wire  measurements  made  with  a Hastings 
Model  HR-2  nondirection al  hot  wire  anemometer. 
The  Hastings  Model  HR-2  is  an  instrument  that 
measures  all  components  of  velocity  in  an  additive 
manner.  When  measurements  of  the  mean  hori- 
zontal flow  were  made  at  a point,  the  values  re- 
ceived on  the  recording  chart  were  not  actually 
true  horizontal  values.  Looking  at  the  flow 
pattern  in  two  dimensions,  the  hot  wire  measures 
the  mean  horizontal  flow  -(vector  A)  and  the 
resolved  component  (vector  D)  of  the  horizontal 
fluctuation  (vector  B)  and  the  vertical  fluctuation 
(vector  C). 


K(zi)=l^j(w'y 


B 
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Table  4. — Wind  'profile  data  jor  wheat  at  Ithaca, 
N.Y.,  July  23,  1960 


Profile 

Height 

Veloc- 

V*lk 

2o 

D 

No. 

ity  1 

Cm. 

Cm.lsec. 

Cm.jsec. 

Cm. 

Cm. 

( 25 

20.  7*  ■ 

50 

34.  6* 

75 

41.  2* 

1 

90 

100 

70 1 
140t 

■123±5 

7.  2 ± . 9 

-77± 

120 

215t 

160 

305t 

240 

380 1 . 

25 

50 

37.  7* 

75 

48.  6* 

2 

90 

100 

120t 

195t 

■131±7 

5.  9±1 

-75±2 

120 

265 1 

160 

345 1 

240 

440 1 

25 

24.  6*  ■ 

50 

52.  4* 

75 

54.  3* 

3 

90 

100 

1401 

215t 

>160  ±6 

7.  4±.  8 

-72±2 

120 

295 1 

160 

400 1 

240 

500 1 . 

25 

25.  6*  ■ 

50 

51.  8* 

75 

64.  5* 

4 

90 

100 

180t 
250 1 

■241  ±3 

14.  2±4 

-60±1 

120 

345t 

160 

470 1 

240 

610t  . 

25 

29.  6* 

50 

89.  9* 

75 

98.  0* 

5 

90 

205t 
280 1 

100 

120 

370t 

160 

490 1 

240 

640 1 

* *Hot  wire  data;  f anemometer  cup  data. 


A typical  hot  wire  trace  may  be  represented 
thus: 


The  mean  velocity  was  determined  by  averaging 
the  readings  at  intervals  of  two-thirds  second  over 
a 30-second  run.  If  the  horizontal  wind,  u,  is 


considered  in  a steady  state  (a  to  b),  the  fluctu- 
ations during  this  period  are  assumed  to  be  the 
fluctuations  of  the  vertical  wind  component,  w' . 
An  averaging  period,  however,  from  0 to  t would 
include  a large  eddy  regime.  Then,  averaging 
over  this  period  would  render  a root  mean  square 
of  the  fluctuations  for  the  horizontal  component, 

V {n'Y,  of  the  large  eddy.  On  the  other  hand,  by 
selecting  for  steady  periods,  a good  estimate  of 

{w'Y  may  be  obtained  {33). 

Use  is  made  of  Taylor’s  Theorem  (54),  which 
employs  a statistical  approach  to  turbulence. 
Using  an  entire  30-second  steady  state  hot  wire 
record  and  displacing  it  two-thirds  of  a second  for 
each  time  increment,  h,  we  assigned  the  whole 
record  equal  weight  in  an  autocorrelation,  Bt, 
where  Rt  is  the  correlation  between  w't^  at  the 
time  ti  and  w«'j4-2/3  at  time  ^1+2/3,  etc.  A good  es- 
timate for  the  mixing  length,  I,  will  then  be: 


where  ^„=time  for  Rt=0. 

A typical  autocorrelation  curve  is  shown  in 
figure  23.  It  is  noticed  that  instead  of  lim  Rt 

t^ca 

— {),Rt  crosses  zero  and  then  oscillates  around  it. 
This  shows  that  the  ideal  case  is  approached  in 
this  evaluation,  but  never  achieved. 

Having  estimated  the  vertical  mixing  length, 
I,  and  the  vertical  component  of  velocity,  we  may 
then  determine  the  eddy  diffusion  coefficient  at 
any  point  in  space  from  equation  121.  Then, 
from  this  relationship  and  the  wind  profile,  we 
may  estimate  the  shearing  stress  (momentum 
exchange)  by: 

T=  pK{Zi)duldz,  or 

CO2  exchange  P=K{Zi)d  C02ldz,  (123) 

where  P is  the  CO2  exchange  rate,  and  d C02ldz 
is  the  gradient  of  CO2. 

Wheat  Investigations 

On  July  23,  the  mast  was  placed  so  that  the  cup 
anemometers  were  at  heights  of  90,  100,  120,  160, 
and  240  cm.  above  the  ground.  The  site  was  the 
same  as  that  used  for  the  previous  measurements 
with  the  wheat  130  cm.  high.  On  this  day,  how- 
ever, a Hastings  Model  HR-2  hot  wire  anemom- 
eter with  an  0.5-second  time  constant  was  used  to 
obtain  mean  wind  velocities  at  heights  of  25,  50, 
and  75  cm.  in  the  plant  canopy. 

Since  there  was  only  one  hot-wire  instrument  to 
scan  the  chosen  individual  levels  in  the  vegetative 
canopy  of  the  wheat,  a “normalizing”  procedure 
had  to  be  adopted  to  permit  a complete  mean 
windspeed  profile  evaluation.  This  consisted  of 
taking  numerous  semisteady  30-second  runs  with 
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Rt 


Z 


(cm.) 


Figuee  23. — Sample  autocorrelation  (Rt)  as^a  function  of 
lag  time  (At)  times  mean  windspeed  («)  for  wheat, 
July  23,  1960. 

the  hot  wire  anemometer  at  a given  level  in  the 
crop  simultaneous  to  numerous  semisteacly  30- 
seconcl  profile  runs  with  the  cup  anemometers 
above  the  crop.  Hot  wire  measurements  at  the 
various  levels  in  the  canopy  then  were  selected 
on  nearly  identical  “above  the  crop  profiles.” 
Five  complete  profiles  were  selected  in  this  way. 

Values  for  the  five  profiles  are  given  in  table  4 
and  figure  24.  There  are  three  observations  of 
immediate  importance:  (1),  the  velocities  at  25 
cm.  are  nearly  the  same  for  all  wind  regimes  cov- 
ered; (2),  the  values  of  dujdz  between  50  and  75 
cm.  are  quite  small  for  each  profile;  and  (3),  the 
point  of  inflection  relative  to  the  crop  surface  ap- 
pears to  be  within  20  cm.  below  the  crop  top. 

From  hot  wire  measurements  made  in  the  wheat 
(130  cm.  high)  and  cup  measurements  made  in  and 
above  the  wheat,  two  profiles  that  will  be  consid- 
ered in  detail  in  this  section  are  plotted  in  figure 
25.  Using  the  method  just  described,  we  calcu- 
late the  values  for  K{zd  for  the  three  heights  of 
25,  50,  and  75  cm.  within  the  crop.  Above  the 
crop,  K{Zi)  was  calculated  from  the  semilog  profile, 
using  equation  120. 

Figure  26  shows  the  distribution  of  K{Zi)  in- 
creasing with  height  above  the  ground  and  in- 
creasing with  windspeed  above  the  crop.  The 
distribution,  K{Zi),  between  75  cm.  and  top  of  the 
crop  was  not  determined,  but  it  is  believed  to  be 


Figure  24. — Wind  profiles  in  wheatfield,  July  23,  1960. 

rather  complex,  based  upon  the  results  found  in 
corn. 

Two  items  of  importance  should  be  pointed  out. 
In  contrast  to  the  findings  reported  by  Penman 
and  Long  {22)  for  wheat,  the  transport  coefficient 
here  is  considerably  larger  than  for  molecular  dif- 
fusion (10^  greater)  and  the  wheat  crop  did  not 
appear  to  “seal  itself”  with  increasing  windspeed; 
i.e.,  the  transport  coefficient  within  the  crop  was 
sensitive  to  windspeed  above  the  crop. 

Calculations  for  shearing  stress  are  given  in 
figure  27.  In  this  case,  shear  does  appear  to  be 
considerably  attenuated  within  the  upper  half  ol 
the  crop. 

Corn  Investigations 

On  August  2,  1960,  wind  profiles  were  taken  in 
Ellis  Hollow  cornfield  {U)  about  30  feet  north- 
west of  the  central  instrument  line.  Northwest 
winds  blew  approximatelv  45°  to  the  corn  rows. 
The  sky  became  cloudy  after  2:30  that  afternoon, 
so  that  near  isothermal  conditions  existed  during 
the  test  period. 

Wind  profiles  were  taken  with  cup  anemometers 
above  the  corn  at  heights  of  250,  300,  and  350  cm. 
The  three  cup  anemometers  were  used  as  a refer- 
ence for  “normalizing”  the  hot  wire  measurements 
of  25,  50,  75,  100,  150,  200,  230,  and  300  cm.  above 
the  ground.  The  crop  averaged  about  240  cin. 
high  and  was  approximately  80  percent  tasseled. 
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Figure  25. — Detailed  wind  profi  les  in  wheatfield,  July  23 
1960. 


At  this  stage  of  development,  there  is  an  unusually 
high  concentration  of  leaves  at  the  top  of  the  plant. 
Subsequent  rapid  vertical  growth  of  the  stem  to 
full  height  eventually  leads  to  a more  even  dis- 
tribution of  leaves  along  the  stems. 

Four  wind  profiles  were  chosen  that  cover  the 
spectrum  of  winds  that  blew  that  afternoon.  Fig- 
ure 28  gives  the  graphical  analyses  of  the  turbulent 
boundary  layer  profiles  above  the  crop.  Too  few 
levels — three — were  taken  to  determine  adequately 
the  profile  parameters,  V*,  Zq,  and  D,  by  any 
statistical  method.  Nonetheless,  by  use  of  judg- 
ment from  past  experience  in  the  graphical  method 
with  mature  corn  {14)  realistic  values  were  arrived 
at  in  this  case.  Pertinent  data  are  given  in  table  5. 

The  measurements  taken  in  and  above  the 
crop  canopy  under  similar  winds  were  plotted 
in  figure  29.  The  solid  lines  were  drawn  in  as  a 
“best  fit”  approximation.  Three  points  should 
be  mentioned  at  this  time: 

(1)  There  is  fairly  good  agreement  between 


z 


Figure  26. — Transfer  coefficient  for  momentum  (Km)  as 
a function  of  height  (z)  in  and  above  a wheat  crop, 
July  23,  1960. 


Table  5. — Wind  'profile  data  taken  above  corn  'with 
cup  anemometers  at  Ithaca,  N.Y.,  August  2,  1960 


Profile 

No. 

Height, 

Zi 

Veloc- 
ity, u 

V*lk 

Zo 

D 

Cm. 

Cm.jsec. 

Cin./sec. 

Cm. 

Cm. 

{ 250 

130 

1 

300 

160 

56.  2 

7.  4 

-180 

[ 350 

180 

f 250 

215 

2 

300 

385 

98.  5 

9.  6 

-130 

[ 350 

350 

250 

360 

3 

300 

405 

144.  8 

12.  8 

-100 

( 350 

435 

250 

420 

4 

300 

470 

177.  7 

17.  0 

-80 

i 350 

505 

the  cup  anemometer  and  hot  wire  at 
the  SOO-cm.  level; 

(2)  duldz  is  maximum  just  below  the  crop 
top;  and 
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z 

(cm.) 


Figure  27. — Drag  (t)  as  a function  of  height  {z)  in  and 
above  a wheat  crop,  July  23,  1960. 

(3)  dujdz  is  quite  small  between  25  and  100 
cm. 

In  explanation  of  the  last  observation,  it  was 
determined  that  the  vertical  component  of  velocity 
was  8 percent  of  the  mean  horizontal  component 
above  the  crop;  however,  this  value  increased 
as  the  height  of  the  measurement  was  decreased. 
At  the  25-cm.  position,  it  was  estimated  that  the 
vertical  component  was  about  30  percent  of  the 
mean  horizontal  flow.  Since  the  vertical  com- 
ponent is  additive,  dujdz  between  25  and  100 
cm.  should  be  greater  than  shown  in  figure  29 
and  table  6.  Therefore,  the  previous  assumption 
of  isotropic  flow  is  less  correct  as  the  measure- 
ments approach  the  ground. 

The  statistical  mixing  lengths  were  calculated 
at  each  point,  using  the  method  described  above. 
The  values  obtained  are  listed  in  table  7. 

The  values  for  above  the  crop  (300  cm.)  are  at 
best  approximate;  however,  they  are  in  the  right 
order  of  magnitude.  Since  the  ratio  of  llz=0A 
within  the  turbulent  boundary  layer  (Karman), 


Z + D (cm.) 


y (cm. /sec.) 


Figure  28. — Log  profiles  of  wind  above  cornfield,  August 
2,  1960. 

Table  6. — Wind  'profile  data  taken  with  hot  wire 
anemo'meter  at  Ithaca,  N.Y.,  Jul'y  2,  1960 


Height, 


Windspeed,  u,  according  to  profile  No. — 


Zi 


(cm.) 


1 

2 

3 

4 

Cm.lsec. 

Cm.lsec. 

Cm.lsec. 

Cm.lsec. 

25 

12.  6 

18.  0 

21.  0 

43.  0 

50 

19.  0 

20.  9 

30.  0 

48.  0 

75 

19.  7 

26.  9 

32.  0 

100 

21.  5 

30.  0 

27.  4 

60.  1 

150 

23.  5 

44.  7 

52.  0 

100.  0 

200 

70.  4 

124.  0 

156.  7 

249.  5 

230 

111.  2 

166.  6 

228.  4 

300.  0 

300 

158.  5 

266.  5 

425.  7 

488.  8 
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Table  7. — Mixing  length  I (z)i,  as  a junction  oj 
height,  at  Ithaca,  N.Y.,  Aug.  2,  1960 


Mixing  length,  I,  for  profile  No.- 


(cm.i 

1 

2 

3 

4 

Cm. 

Cm. 

Cm. 

Cm. 

25 

21.  0 

25.  0 

25.  5 

30.  0 

50 

31.  0 

33.  0 

36.  0 

47.  0 

75 

39.  5 

45.  0 

47.  5 

100 

50.  5 

50.  0 

53.  0 

55.  5 

150 

44.  0 

53.  0 

78.  0 

93.  0 

200 

75.  5 

75.  5 

100.  0 

no.  0 

230 

183.  0 

196.  0 

205.  0 

220.  0 

300 

63.  0 

74.  5 

140.  0 

170.  0 

one  would  expect  close  to  the  same  value  at  300 
cm.  Using  the  ground  as  the  reference  plane,  we 
obtain  Ijz  range  from  0.21  to  0.565.  With  the 
zero-plane  displacement  given  as  the  reference, 
the  values  range  from  0.485  to  0.71. 

The  mixing  length  values  at  25  cm.  are  con- 
sidered to  be  too  large.  This  is  attributed  to  the 
high  fraction  of  vertical  component  of  velocity 
when  calculating  the  mean  horizontal  flow. 
Within  this  region,  Ijz  is  approximately  1. 

The  values  for  Ijz  are  fairly  constant  between 
50  and  200  cm.  They  are  in  the  magnitude  of  0.5. 
However,  the  most  mysterious  values  were  ob- 
tained at  a point  10  cm.  below  the  crop  top 
(2=230).  At  this  point,  the  mixing  lengths  are 
larger  than  those  values  obtained  in  the  boundary 
layer  above  the  crop,  and  are  therefore  suspect. 
The  observation  made  while  the  data  were  being 
analyzed  indicates  that  the  reason  for  these  high 
values  originates  from  the  abnormally  long  period 
of  the  cycle  about  the  mean  velocity. 

At  the  position  of  230  cm.,  the  frequency  of  one 
cycle  about  the  mean  was  from  three  to  five  times 
the  ones  previously  encountered.  It  is  suggested 
that  this  may  be  due  to  the  elastic  waving  of  the 
stalk  and  the  upper  leaves  superimposed  upon  the 
normal  frequency.  The  longer  cycles  were  respon- 
sible for  a much  less  rapid  drop  in  the  auto- 
correlation with  time,  and  hence  a larger  value  for: 


Values  for  the  eddy  diffusion  coefficient,  K{Zi), 
were  calculated  from  equation  121  (table  8).  They 
were  then  plotted  in  figure  30.  The  dashed  lines 
in  figure  30  represent  Kizj  as  calculated  from  the 
semilog  profile,  and  the  solid  lines  are  from  the 
hot  wire  data.  The  shearing  stress  was  calculated 
from  equation  122  (table  9).  Figure  31  shows  the 
distribution  of  the  shearing  stress. 


Figure  29. — Wind  profiles  in  and  above  cornfield,  August 
2,  1960. 

Table  S.— Transfer  coefficient,  K{zi),  as  a function 
oj  height  at  Ithaca,  N.Y.,  August  2,  1960 


Height, 

Zi  (cm.) 

K{zi)  for  profile  No. — 

1 

2 

3 

4 

Cm.^/sec. 

Cm.^lsec. 

Cm. sec. 

Cm.^lsec. 

25 

94 

208 

298 

882 

50 

190 

287 

428 

1,  349 

75 

181 

486 

561 

100 

359 

540 

1,  039 

1,  576 

150 

339 

1,  211 

2,  223 

2,  818 

200 

1,  306 

1,  910 

2,  700 

3,  894 

230 

5,  106 

5,  606 

6,  724 

7,  854 

300 

1,  796 

2,  712 

5,  250 

7,  208 

SUMMARY 

The  determination  of  Km  in  and  above  the 
momentum  sink  elements  using  hot  wire  anemom- 
eter data  affords  an  independent  check  on  the 
semilog  profile  method  for  determining  Km,  and 
also  permits  a distribution  analysis  of  where  the 
momentum  is  being  transferred  in  the  vegetative 
sink.  Of  particular  interest  here  is  the  com- 
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Z 


(cm.) 


Figuke  30. — Distribution  of  transfer  coefficient  for  momentum  {Kvi)  in  and  above  cornfield,  August  2,  1960. 


Table  9. — Shearing  stress  as  a function  of  height 
at  Ithaca,  N.Y.,  August  2,  1960 


Height, 

Zi  (cm.) 

Shearing  stress  for  profile  N 

0. — 

1 

2 

3 

4 

Dynesjcm.^ 

Dynes/cm.^ 

Dyneslcm.^ 

Dynesjcm.^ 

25 

0.  038 

0.  083 

0.  119 

0.  353 

50 

. 076 

. 115 

. 171 

. 540 

75 

. 073 

. 194 

. 224 

100 

. 143 

. 227 

. 499 

. 946 

150 

. 203 

. 872 

2.  13 

3.  38 

200 

1.  88 

3.  21 

5.  51 

9.  35 

230 

9.  19 

13.  45 

20.  17 

28.  27 

300 

1.  08 

2.  28 

5.  04 

8.  65 

parison  of  the  Km  values  derived  by  the  two 
methods  in  the  turbulent  boundary  layer  above 
the  corn  (Zi=S00  cm.)  in  figure  30.  In  three  out 
of  the  four  profiles  the  hot  wire  values  are  higher 
but  of  the  same  magnitude  as  the  log  profile 
values.  The  largest  divergence  occurs  at  the 
lowest  windspeed,  although  there  appears  to  be 
no  consistent  relation  between  divergence  and 
windspeed.  If  the  largest  divergence  at  slow 
mndspeed  is  indeed  real,  it  should  be  anticipated 
from  the  following  factors : 

(1)  Rotating  cup  anemometers  are  subject 
to  greater  errors  at  slow  windspeeds  ; 

(2)  There  are  errors  in  estimating  'firee 
convection”  from  the  log  profile  caused 
by  thermal  gradients ; the  latter  are  more 
pronounced  at  low  windspeeds. 
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It  should  be  pointed  out,  however,  that  tem- 
perature gradients  were  extremely  small  in  the 
turbulent  boundary  layer  over  the  corn  (tempera- 
ture gradients  may  have  been  significant  within 
the  canopy,  however),  and  that  the  slowest  wind 
measured  above  the  crop  was  well  above  the 
stalling  speed  of  the  cup  anemometers. 

The  distribution  of  shear  in  figure  31  can  be 
considered.  Realizing  that  the  techniques  used 
and  the  assumptions  made  are  open  to  criticism, 
and  that  the  results  should  be  viewed  as  only 


semiquantitative,  we  should  point  out  that  the 
origin  of  shear  from  the  log  profiles  falls  on 
‘‘planes”  between  80  and  180  cm.  above  the 
ground  {—D  values.)  On  the  other  hand,  the 
hot  wire  evaluation  of  shear  within  the  canopy 
demonstrates  a “volume”  distribution  of  shear 
in  a complex  but  reasonable  fashion.  The  hot 
wire  quantitative  values  are  obviously  too  high 
if  the  assumption  is  made  that  the  total  shear  of 
the  “surface”  is  accurately  evaluated  by  the  log 
wind  profile  characteristics  above  the  “surface.” 


Z (cm.) 


Figure  31. — Distribution  of  drag  (t)  in  and  above  cornfield,  August  2,  1960. 


SUMMARY  AND  CONCLUSIONS 


1.  Two  theoretical  models  were  developed  for 
wind  flow  characteristics  in  plant  canopies  based 
upon  a general  type  differential  equation  for  the 
vdnd  profile  both  inside  and  above  the  vegetative 
canopy  at  the  earth’s  surface. 

2.  A testing  of  the  theoretical  models  depended 
upon  a breakdown  of  the  wind  structure  into 
quasi-steady  and  transient  states.  It  was  assumed 
that  the  small  scale  eddies  caused  by  the  micro- 
structure of  the  surfaces  (plant  surfaces)  are 
responsible  for  vertical  transport.  For  a begin- 
ning, the  simpler  case  was  studied  experimentally. 
Wind  flow  in  and  above  the  vegetation  at  the 
quasi-steady  state  produced  trends  predicted  by 
the  turbulent  flow  model. 

3.  The  turbulent  canopy  flow  model  postulated 
a linear  growth  of  diffusivity  within  the  turbulent 
canopy  flow  layer  rather  than  the  log  velocity 
profile.  This  led  to  a realistic  curvature  reversal 
of  the  velocity  profile  deep  in  the  canopy  flow 


layer.  This  must  be  a basic  characteristic  of  the 
turbulent  wall  flow^  within  the  roughness  elements. 

4.  Studies  of  the  log  velocity  wind  profile 
characteristic  over  various  kinds  of  vegetation 
revealed  a complex  coupling  betw^een  windspeed, 
roughness  length,  friction  velocity,  and  effective 
displacement  and  the  surface  geometric  and 
elastic  properties.  In  twm  cases,  the  friction 
velocity  remained  constant  and  independent  of 
windspeed.  To  facilitate  these  and  future  studies, 
a machine  computation  method  was  developed 
that  includes  an  evaluation  of  the  standard  errors 
of  the  logarithmic  wind  profile  parameters. 

5.  An  evaluation  of  the  transfer  coefficient  and 
shear  within  the  plant  canopy  revealed  a reason- 
able picture  of  their  spatial  distribution  within 
the  “crop  volume.”  This  serves  to  emphasize 
that  the  concept  of  the  roughness  length,  although 
very  useful,  is  in  reality  an  artifact. 
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